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SUMMARY 


This report presents the derivation of a frequency domain theory and working equations for 
radiation of propeller harmonic noise in the presence of angular inflow The 
the acoustic ^analogy generalized for a moving medium by M. E. 

In applying the acoustic analogy, integration over the tangential coordinate of the source regio 
is performed numerically, permitting the equations to be solved without approximation for ^y 
decree of angular inflow. All source terms are included. The thickness source is represent^ 
by a monopole distribution that is unsteady by virtue of varying blade s^tion re ative sp^. 
Loading, represented by dipole distributions, may be unsteady and includes radial, tangenti^, 
and axial components. The quadrupole term includes all nine components and ^so may be 
unsteady. General equations are written for sources on the actual blade surfaces; these give t e 
sound field exactly, assuming that the source terms are known and that there are no surfaces 
present other than the blades. However, for most purposes, a mean surface approximation will 
be adequate. In this case, the sources may be placed on the camber surface (or the chord line). 
For far field calculations, special forms of the equations are given to reduce computer running 
time. Angle of the propeller inflow is specified in terms of yaw, pitch, and roll angles of the 
aircraft. Since these can be arbitrarily large, the analysis applies with equal accuracy to 
propellers and helicopter rotors. 

The quadrupole derivation has been carried only far enough to show feasibility of the 
numerical approach. However, for thickness and loading, the derivation is given in complete 
detail with working equations for near and far field. Explicit formulas are presented for 
computation of source elements, evaluation of Green’s functions, and location of observer points 
in various visual and retarded coordinate systems. The resulting computer program, called 
WOBBLE, has been written in FORTRAN and follows the notation of this report very closely. 
Only about 60 lines of code are required to compute the linear radiation terms exactly. 

The new theory was explored to establish the effects of varying inflow angle on axial and 
circumferential directivity. Also, parametric studies were performed to evaluate various 
phenomena outside the capabilities of earlier theories, such as an unsteady thickness e ec . 
Widity of the theory was established by comparison with test data from conventional propellers 
and Prop-Fans under a variety of operating conditions and inflow angles. Agreement with the 
propeller data is excellent and provides a satisfying verification of the theory. Agreement with 
the Prop-Fan data is not as good. Since the acoustic equations are exact, the difficulty must be 
in estimating the aerodynamic blade loading and not with the acoustic radiation equations For 
Prop-Fans, tip flow is influenced by vortex loading which can dominate the noise at subsonic 
tip speeds but is of secondary importance for aerodynamic performance. 


SECTION 1 
INTRODUCTION 


It has been known for some time that propeller noise can be strongly affected by any inflow 
that l^ds to unsteady blade loading. Experimental and theoretical resear ch in„this area 
intensified in the 1980 s primarily as a result of the advanced turboprop program at NASA. It 
was believed that noise due to angular inflow could be predicted accurately by computing the 
blade loading from an unsteady aerodynamic analysis and then using an axial flow model for the 
noise radiation. However, in 1990 Mani (Ref. I) showed via a first order analysis S 
flow Mach number could have a more important impact on the noise Mn wouldlhe un 
loading. Krejsa (Ref. 2) extended Mani’s work to all orders and showed an impressive 
agreement with test results. The work at Hamilton Standard grew out of discussions with the 
noise group at NASA Lewis and with Mani. The first contribution from Hamilton Standard 
was a paper by Hanson (Ref. 3) which showed that noise in the far field coulJ be described in 
terms of "wobbling modes" whose simple form clarified the physics of the noise generation 
process. The noise increase above the levels predicted with unsteady loading but axial inflow 
was traced directly to the increase in the Mach number of the blade sources with respect to the 
observer. The second contribution from Hamilton Standard is the work reported herein. 

The objective of this work was to develop the noise theory for the near field via a direct 
frequency domain analysis valid when the propeller shaft is inclined at ainy angle to the flight 
direction, as suggested by Figure 1. Geometrical approximations of the previous freque^^^^ 
domain methods have been eliminated so that the radiation equations are exact for specified 
linear sources. This report presents the derivation of the theory, results of exercising the theory 
in parametric studies, and a validation of the theory by comparison with test data. Note that the 
analysis applies strictly to pure ongular inflow, not to distorted inflow. 

TJnfike previous works by Hanson which place sources on the propeller advance surface, the 
new equations place sources on the blade surfaces assuming this source information (blade 
pressure distribution and geometry) is available from a finite element aerodynamic solution 
However, guidance is also given for use of the thickness and loading sources specified in terms 
of blade section thickness and lift distributions placed on a mean camber or chordal surface. In 
the near field formulation, loading variation may be represented directly as a time history if 
desired, eliminating summations over spinning modes. The far field formulation is an exLct 
extension of the near field model and uses identical source input. Tangential integration is 
performed analytically, leading to Bessel functions as in earlier far field formulations and 
reproducing the "wobbling mode" behavior first revealed in Ref. 3. However, the chordwise 
integration is left in numerical form to retain the exact nature of the solution. Although the 
derivation is somewhat complex, the working equations are simple and easy to interpret. 

The FORTRAN source listing, compiled code for the IBM PC, and associated input and 
output data sets are available from the Contract Manager, Mr. Bruce Clark, at the NASA Lewis 
Research Center, 21000 Brookpark Road, Cleveland, Ohio 44135. Telephone number 216-433- 
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SECTION 2 
BACKGROUND 


For simplicity, the earliest analyses of propeller noise (in the 1930’s and 40’s) moiled the 
blades as body forces in an infinite medium with no other surfaces and no mean flow. Forward 
fl^^ht e^were added by Garrick and Watkins (ref. 4) in the 1950’s^ 
propellers in an ideal environment, predictions from these theones compared favorably with test 
data However, propellers of practical interest are installed on aircraft with reflecting surfaces 
and with inflow that generates unsteady loading. Under these conditions, agreement between 
noise theory and test data was often disappointing. Unsteady blade loading caused by non-ideal 
inflow was believed to be the problem. Progress toward representing the true flow environment 
was slow. Noise analyses for unsteady blade loading without forward flight were studied 
extensively in the 60’s, most notably by Morse and Ingard ^ef. 5), Lowson ^ef. 6), and 
Wright (Ref. 7). Forward flight was added to these models in the 70’ s and 80 s by Hanson 

(Refs. 8, 9). 


In the 80’s Farassat refined the modeling in his time domain methods for propellers and 
helicopter rotors (WOPWOP, for example) by accounting for a general motion of the rotor 
including that where the rotor axis was tilted with respect to the flight direction This ^pect 
of his method was not explored extensively but it is documented in a paper by Padula and Block 
(Ref 10) It wasn’t until the late 1980’s that the angular inflow problem was given senous 
consideration. In 1988, R. Stuff published the first frequency domain theory for rotating point 
loading sources with inflow at an angle to the axis of rotation (Ref. 11). Soon thereafter, work 
on angular inflow effects was sponsored by NASA-Lewis with Mani (Ref. 1) and Krej^ (Re . 
2) as principal investigators. These and Hanson’s wobbling mode theory (Ref. 3) were fw field 
methods that revealed the importance of the cross flow Mach number but still included some 
geometrical approximations. In the far field methods, the intention over the source in the 
tangential direction was done analytically, leading to Bessel functions in the radiation equations. 
Meanwhile, Envia had been working on a near field method (Ref 12) in which the 
circumferential integration is done analytically as in the earlier Garrick and Watkins work (Ref 
4). In order to explain near field observations on the Prop-Fan Test Assessment aircraft, 
Hanson and Farassat both generated near field methods capable of handling angular m^w. 
Hanson’s method (Ref. 13) was a simplified version of the present theory and Farassat s (Ref 
14) was a program in which he reinstated angular inflow capability in his time domain code. 


The present work returns to that same idea for the tangential integration but with the precise 
treatment of blade geometry and inflow angle that was previously available only with the direct 
time domain methods of Farassat and others. 
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SECTION 3 
THEORY DERIVATION 


The problem addressed in this section is noise of a propeller operating in a uniformly 
flowing, unbounded medium containing no surfaces or bodies other than the propeller. Although 
this is treated with complete generality, there are limitations that should be explained. In 
particular, note that angular inflow does not imply distorted inflow. Thus, for example, noise 
caused by wakes of an upstream wing could be calculated but the effect of refraction of the 
sound in the wake shear flow and reflection from the wing would not be represented. This 
section starts with a review of Goldstein’s version of the acoustic analogy (Ref. 15), which is 
the starting point for development of the propeller noise theory. First, general equations are 
derived in the form of continuous integrals over the source. These are then discretized in a form 
directly applicable to computer coding with finite elements of blade area. Finally, special forms 
for a faster running code are derived for the far field. Some of the details on coordinates. 
Green’s functions, and source specification are relegated to Appendices A, B, C and E. 
Appendix D gives a detailed treatment of observer coordinate systems to facilitate prediction at 
points specifed in visual or retarded coordinates and in systems aligned with the flight direction 
or aligned with the propeller axis. Finally, Appendix F provides a list of symbols. 


THEORETICAL FOUNDATION - GOLDSTEIN'S ACOUSTIC ANALOGY 

Goldstein extended Ffowcs-Williams & Hawkings acoustic analogy formulation (with moving 
surfaces) to the moving medium case. For reference purposes, Goldstein’s formulas are 
presented below with slight changes in notation for the current application. His Equation 13 
gives the acoustic density disturbance p ' at the field point x as the sum of 3 integrals: 

clAx,)--lll^^pV„^clA(y)dr 


+ 




[M f^^dA(y)dr 
J -T J A(r) dy. 

f^J^dydr 

J -rj V(r) dy.dyj 


( 1 ) 


where c„ is the ambient speed of sound. At sufficient distance from the source, cjp ' is equal 
to the acoustic pressure p '. Source coordinates are y=(y„y 2 ,yj) and source time is r. In each 
of the 3 terms, the outer integral is over a range of source time -T<t<T large enough to 
include all contributions to the noise signal at observer time t. 


The first term in Equation 1 is the thickness (or volume displacement) noise in which is 
the ambient density and is the velocity of the surface normal to itself relative to the fluid. 
Interpretation of the thickness source is given in Figure 2 which shows that Vf^dAdr is the 
volume displaced by the surface element dA in the time increment dr. (Note: in computing 
for Equation 1, the surface normal vector points outside the blade, opposite Goldstein’s 
convention.) G is a Green’s function and D/Dt is the convective derivative. The inner integral 
in Equation 1 is over the surface area which, as indicated by the notation, may be moving. 

The second term in Equation 1 is the loading noise in which / are the vector components 
of the force per unit area on the surface. Thus, f^dA is the element of force on the area element 
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dA. In this term, summation of i over the 3 space indices is implied via the Einstein 
summation convention. 

The third line in Equation 1 is the quadrupole term, which accounts for all non-linear effects 
as an integral over the source volume V. is Ligh thill’s stress tensor 

7^. = p V,. V. + 5 .. [(p -pj - cIp^ - e,j (2) 


in which p is the density, v, are the components of velocity disturbance, p-p„ = p' is the 
pressure disturbance, and 5,y is the Kroneker delta (1 for i=j, 0 for i ^j). is a viscous term 
defined in Ref. 15. To evaluate 7^ as a source for noise calculations, the first 2 terms in 
Equation 2 could be obtained from an Euler solution. The viscous effect in the last term could 
be obtained from a Navier-Stokes solution or from a boundary layer code. 

A quote from Goldstein (Ref. 15) indicates the significance of Equation 1: 

"This is the fundamental equation governing the generation of sound in the presence of 
solid boundaries. It is an exact expression and (provided the integrals converge) applies 
to any region which is bounded by impermeable surfaces in arbitrary motion. In the 
acoustic analogy approach we assume that the stress tensor Ty' and the applied force f 
can be modeled mathematically or determined experimentally. Then, the right-hand side 
of this equation is known and the density fluctuations in the sound field can be 
calculated." 

In the following subsection, these integrals will be applied to the propeller problem without 
approximation. Then discrete forms of the integrals and far field equations are developed. A 
mean surface representation of the blade geometry is appropriate for most engineering work and 
is described in Appendix C. 


DERIVATION OF GENERAL NOISE RADIATION EQUATIONS 

The key to developing efficient forms of Equation 1 for propeller noise calculation is 
choosing the best coordinate systems for the problem. This section therefore begins with a 
discussion of the propeller coordinates and how they appear in the Green’s function. Then 
radiation integrals are developed for each of the 3 sources. 

Coordinate Systems and the Green’s Function 

This is a moving medium analysis so that the aircraft is fixed with respect to the coordinate 
system and the undisturbed fluid moves past it at a uniform Mach number, M. Observer, or 
field point, coordinates are denoted by x = (x,y,z), a right hand system with the x-axis pointing 
in the flight direction. The y-axis could lie in any direction but usually will be horizontal, 
pointing to the right of the flight path. This system, sketched in Figure 1, conforms with 
standard flight dynamics conventions (Ref. 16) with the z-axis pointing down. Various source 
coordinates, denoted in general form as y=(yi,y 2 >yi)> will be used but the basic set, 
y=(x„,y„,zj, is coincident with the flight coordinates as shown in Figure 1. 
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Given this system, the free space, moving medium, time dependent Green’s function given, 
for example, by Blokhintsev (Ref. 17) is 


G = 


8(t - T -ale) 
4x5 


(3) 


where 5 is called the amplitude radius 

^ = sjix-xf +^\(y-yf+{z-zf] 
with and a is the phase radius 

The specific form of the convective derivative should also be noted: 

D d ,, a 

Dt 3t " dx^ 

In these expressions for the phase radius and the convective derivative, the sign of the Mach 
number M corresponds to the present use of j:-axes pointing upstream. 

For the case where the propeller is tilted with respect to the flight direction, it is much more 
convenient for the source to use cylindrical coordinates aligned with the shaft because then the 
propeller rotation involves only one of the coordinates. Source coordinates, defined in Figure 
3, can be obtained by rotations from the reference set through the yaw, pitch, and roll angles. 
The transformations are derived in Appendix A. However, it is sufficient for present purposes 
to note that x„,y„,Zo can be expressed as functions of r„,<j>o,K so that all of the derivatives of 
the Green’s function required by Equation 1 can be performed analytically. 

Thickness Noise 

For simplicity, we first analyze the noise from a rotor with only one blade; at the end of this 
subsection the other blades are accounted for by superposition. The pressure disturbance can 
be written as the sum of the 3 terms in Equation 1: 

cl p' =p'{x,t)= p^x,t) +p^(x,0 +PQ{x,t) ('7) 

so that the thickness noise is 

In this expression, Vf^dA is the volume displaced per unit time by the area element dA due 
to the motion of the blade. If this is considered a source element, then the double integral 
simply represents the summation over the surface and over source time of these elements 
weighted by the Green’s function derivative. The notation A(t) indicates that the blades are 
moving. We prefer to do the surface integral in blade-fixed coordinates, which are defined in 
Figure 4. It can be seen that blade sections are defined by "cutting" the blade with cylindrical 



(4) 

(5) 
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surfaces of radius r„=r,. The figure shows a typical section "unwrapped" onto a plane in which 
the coordinates are and fs<t>s' this system, motion of the blade involves only the (j)^ 
coordinate. Thus, the area element is located at r„, <f)^, xj = r„ 0,+Or, x,. Now, we can 
rewrite Equation 8 as 


p^(x,t) = 


V. 


N 


DG 

~Dt 


dA(y)dT 




( 9 ) 


with the new notation on the brackets indicating that the Green’s function derivative is to be 
evaluated analytically in terms of and evaluated at the surface values shown. may 

be regarded as a function of 4>^, and x„' only, in which case the notation in Equation 9 will 
provide the correct time dependence. With this notation, the A on the inner integral is no 
longer dependent on source time. Since the integration limits do not depend on the variables, 
the order of integration may be switched as follows 


PAx,t) = ~Po\^\_ 


DG 

Dt 


drdAiy) 




( 10 ) 


where r„=r^ and x^^x^ have been dropped from the notation and are implicit henceforth. 
Since pAx.t) is periodic at frequency WIt, it can be represented by the Fourier series 


n=-oo 


( 11 ) 


in which the Fourier coefficients would ordinarily be evaluated according to 

J neriod 


( 12 ) 


However, when Equation 12 is applied to Equation 10, indeterminate limits (t = -7 to 7) remain 
on the source time integral. To avoid this, we adopt the following strategy: consider one 
revolution of the rotor, e.g. for 0 < t < 2-ir/Q, because we know that this will produce exactly one 
period of the noise signal, say from time t, to r,+2ir/fl. (We don’t need to know tj.) The 
resulting noise pulse is 

prixd) =PT(x,t), ( 13 ) 

= 0 , elsewhere 

Now, with the restricted range of the t integral, the Fourier coefficient can be computed from 
a modified form of Equation 12: 
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( 14 ) 


Prni^)=-^r Prix,t)e'''^'dt 
Itt } -« 


To apply this, we insert the time dependent Green’s function from Equation 3 into Equation 10 
to obtain the following 


Prni^) 




27 T/n 


P yKj — 

r' o N 


D 


— f ”6(r-r-a/cje''-"' 
47 t5J-» 


dt 


drdA (15) 


6^4>*Qt 


which, upon performing the /-integration and denoting the harmonic Green’s function by 


G. = 


4tt5 


(16) 


becomes 


P (X) = -^ f f \p^V^—(G„e‘’‘^j 

27t J ^Jo " ' 


drdA 




(17) 


Application of Equation 6 leads to 




r r 2WQ 


u. 

P 

^ o N 


_ 


dG 

inQG^-cM — - 
dx. 


drdA 


(18) 


+Or 


We now change the variable of integration on the inner integral using to get 


rnv 


f ± 

[2. 

p 

dG 

c M " - inSlG 


J 2'3T . 

1 0 

^ o N 

" dx 

o 



e'^^-dA 


(19) 




Because <^>o=</>o^ we can now simply drop the prime on 4>o and drop the notation outside the 
square brackets. Recall that the development so far has dealt with only one blade. For 
interaction with angular inflow or fixed flow distortion, each blade experiences identical loading 
and therefore produces identical noise signals, except for time lags according to blade position 
in the hub. For B equally spaced blades, only the blade passing frequency harmonics (n—mB) 
survive and these are multiplied by B. With this change, the thickness noise integral is 
established: 




C^M. 


dG. 




JmB6 


d4>M 


( 20 ) 
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For programming purposes, a non-dimensional form of Equation 20 is preferred. If we non- 
dimensionalize lengths by propeller tip radius r,, area by time by fl, and pressure by the 
ambient pressure (which equals where 7 is the ratio of specific heats, 1.4 for air), t e 

result is 




dG 
M — ^ 


dA 


( 21 ) 


where M„=V^c„, k„=mBM„ and M,=Ur/c^ (the tip rotational Mach number). The mB 
subscripts have been simplifi^ to m so that the Green’s function is now 

G =— 

" 4x5 


Note that this normalization has changed the definitions of some of the symbols but the meanings 
should be obvious from the context. For example, dA is now the elemental area divided by r, . 

Equation 21 is an exact expression for the thickness noise and is straightforward to evaluate, 
given say the output of an Euler program for propeller performance. Such a program could 
provide the actual <i 4 ’s and the associated with each dA as defined in Figure 4. 

As shown in Appendix A, the Green’s function derivatives can be performed analytically from 
Equations 3-5 and the appropriate coordinate transformations. Then they and the Green s 
function itself can be evaluat^ numerically at r^, If the area elements are counted by 

two indices, a triple loop is required to evaluate Equation 21. 


Loading Noise 

Since this term involves a surface integral, the scheme above for thickness noise can be 
followed and the result can be found quickly. From Equations 1 and 7 the loading noise is 

Time dependence of the blade loading will be expressed non-dimensionally as / (Ut). An 
impiortant distinction regarding the Green’s function derivatives in the thickness and loading 
integrals should be made at this point. For thickness, the derivative dG/dx„ is in the flight 
direction, making it especially easy to evaluate. For loading noise in Equation 23, we should 
first regard all the components of ^ dG/dy^ together. Since the repieated indices imply a sum, 
this is /• VG in vector notation. Because this is a scalar, it can be evaluated in any coordinate 
system. The obvious choice is to use propeller coordinates since they are most convenient for 
representation of loading and since the Green’s function derivatives (one of which is dG/dx ^' ) 
can be performed analytically given Equations 3-5 and the equations for transformation of 

coordinates. 
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By the same arguments used to arrive at Equation 10 from Equation 8, we find 




^yi 


drdAiy) 


( 24 ) 




Then, by using the strategy of selecting one rotation in source time from 0 to 2x/Q, computing 
the loading harmonic from the formula analogous to Equation 14, and inserting Equation 3 for 
the Green’s function, we find 


n 




2t/Q 


^yi 




drdA ( 25 ) 


0^=0^ +Or 


By using Equation 16 for the harmonic Green’s function and accounting for multiple blades this 
leads to 




BQ 

2t 


u 


2ir/n 


//Qt) 

^y. 


e'”'^^drdA 




(26) 


By substituting <f)J=^^+UT as in the thickness noise derivation, we find 



-imB4> 




dG. 


^yt 


mB ji/ia 

e d(l>^ oA 


(27) 


Comments similar to those given in conjunction with Equation 21 apply. Any aerodynamic 
Performance code should give elements of blade loading, fidA , needed for Equation 27. The 
load elements are to be defined in cylindrical coordinates: The Green’s 

function derivatives are to be performed analytically based on the geometric properties of the 
coordinate rotations and then evaluated numerically at the source position (f>^, x/. To be 
more specific about the derivatives, d/dy,=d/dr„, d/dy 2 =(l/rjd/d<l>„, d/dyj=d/dx}. ' 

The non-dimensional form of Equation 27 is found with the same normalization used for 
Equation 21 and with / and normalized by ambient pressure to obtain 





(28) 


As with the thickness noise, the B has been dropped from the subscripts on P and G for 
simplicity. 

Quadrupole Noise 

Because this term requires a volume integral, the derivation proceeds somewhat differently. 
Furthermore, there is no way to superpose individual blade fields so we must deal with the entire 
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propeller field at once. 


From Equations 1 and 7 the quadrupole integral is 



d^G 

dy^dyj 


dydr 


( 29 ) 


Undoubtedly, we will want to evaluate the integrand in propeller coordinates (rather *an flight 

coordinates) so that we can take the y,’s to be and write Ty' -Ty ^ 

the quadrupole source will probably come from an aerodynamic performance code with field 
quantities in blade-fixed coordinates. Say that in the blade system we are given y - 
f r) and let this be defined everywhere (including inside the blade where it is defined 

to** be 0)! In these terms, the desired source for Equation 29 is 


7^(r„ ,r) = f;.(r„, <l>„ - QtX, r) 


(30) 


since <^o ~ 3" fl'r* Now the quadrupole integral is 

T r ^ ^ ^ d^G 


Po(X,t)= f U all 

J * space 


dyM 


dydr 


(31) 


Since the limits on the inner integral now don’t depend on t, the order of integration c^ be 
changed. As with the thickness integral, we pick only one revolution in source time to obtain 
one cycle of the signal in observer time 


r° Ur,.<>.-ary„T) 

* space 


d^G 

3>'.3y, 


drdy 


(32) 


Then, the harmonic formula analogous to Equation 14 for the quadrupole gives 



e 

dyidyj 


dr rJ4>^dXadr^ 


(33) 


where the spatial integration variables are now shown explicitly. 

Equation 33 provides the means to evaluate quadrupole noise h^monics given the flow field 
information from a CFD code. There are various ways to deal with the quadruple integral but 
further work will be postponed until the actual form of the flow field data is specified. In 
particular, note that the formula does not yet recognize that only the noise harmonics n=mB 
will be non-zero. This can be remedied by recognizing that the flow in adjacent blade gaps is 
the same but lagged in time by lir/UB. When this is accounted for, the range of either the t 
integral or the <^>o integral can be reduced by the factor 1/5. 

The formulas derived above for the thickness, loading, and quadrupole noise components are 
expressed in the classic notation of continuum mechanics and were denved without geometncal 
approximations. Practical applications will involve discretization of the area elements. This is 
described in the following subsection. Other approximations such as placement of sources on 
mean surfaces (rather than on the upper and lower airfoil surfaces) and far field forms may also 
be desirable. These are discussed later in this subsection and Appendix C. 
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DISCRETIZATION OF AREA ELEMENTS 
and 

EQUATIONS FOR NEAR FIELD CALCULATIONS 


Equations 21 and 28 from the previous section represent thickness and loading noise 
harmonics as continuous integrals over the blade surface area. The program to compute these 
integrals breaks the area into discrete elements and steps through these using counters for radius 
and chordwise distance. This section provides the notation for these discrete forms and makes 
a few comments regarding their application. The intent was to generate one subroutine for near 
field theory that can accept input data from a full surface source description (upper and lower 
airfoil surfaces) or from single surface representations such as the camber surface, the chordal 
surface, or even the advance surface as in the earlier Hamilton Standard methods. The key to 
this generality is identifying the fundamental source elements that could be generated by any 
system. 


Radial and chordwise counters will be fi and v, respectively. Thus, discretization of the 
area integrals is denoted by \ ^dA-^ ,. For thickness noise, the fundamental source element 

is Mf/dA. In discrete form this thickness source element will be defined as , where 

(34) 


where, as defined in the previous section, M^=V^/c„ and is the velocity of the area element 
normal to itself relative to the fluid. (Various means for evaluation of are given in 

Appendix B.) The fundamental loading source element from Equation 28 is f^dA whose 
discrete form is defined to be 






(35) 


Both and , are nondimensional because the area elements are normalized by and 
fi is force per unit area on the blade surface normalized by the ambient pressure. The i index 
counting the load components is interpreted as 


F, 

= F 

l/>.r 


F, 

= F. 

2ii,r 



= F, 


X 


(36) 


With these definitions for the source elements, the discrete area forms for thickness and loading 
noise harmonics can be written down immediately: 


/■..W = tBE 
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(37) 


and 
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Except for the converting the 4>o integrals to discrete form, these are the working forms for the 
near field integrals. For far field calculations, these integrals are performed analytically, leading 
to Bessel functions as with previous Hamilton Standard formulations. However, unlike the 
previous formulations, no approximations regarding geometry are necessary. 

Working forms for the Green’s function and its derivatives are given in Appendix A. 
Information on computing the source terms , and F),, , is given in Appiendix B. Blade 
geometry enters through location of the area elements. These are specified at radii = r^, 
which is an input array for the computer program. At each radius, the axial and tangential 
positions are specified by 2 dimensional arrays x^' = , and <f>^ = and x^^ , are used 

in computing the Green’s function and derivatives; is used directly in ^nations 37 and 38. 
This blade geometry information could be obtained from a finite element aerodynamic 
calculation, if the full surface coordinates are to be used, or from the mean surface 
approximation given in Appendix C. 


DERIVATION OF FAR FIELD RADIATION FORMULAS 

For observers far from the propeller, the far field version of the Green’s function derived in 
Appendix E can be used to generate a faster computational algorithm. Furthermore, the 
formulas eliminate numerical precision problems that could arise from subtracting large, nearly 
equal numbers when integrating over source coordinates. This section gives the derivation of 
far field integrals for thickness and loading noise based on Equations 37 and 38 above. 
Observer coordinates are expressed using notation introduced in Appendix D. 


Thickness Noise 

As shown in Appendix E, the term in brackets in Equation 37 can be written 


dG 

M-^ -ikG = 
dx„ " 


T ■ 

Uf 


(39) 


and the far field Green’s function can be written 


Jk (r-Sx'.) 

” 4x5. 


(40) 


where the upper/lower sign applies to right/left hand rotation of the propeller and where dj is 
the Doppler factor: 

J^=l-Mcos0 (41) 


Here, the observer position is defined in terms of the variables r, S„, S^, and S^, which were 
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defined in Appendix D. These variables can be computed either from visual or retarded 
coordinates; however, the radiation integrals are more easily interpreted in physical terms using 
the latter, as was shown in Ref. 3. 

With the appropriate substitutions, the far field version of Equation 37 becomes 


■ 4x5,4 - 




(42) 


where 


T,,’ 2xJo 


(43) 


The remaining work is to evaluate the tangential integral „ which we do by bringing out the 
explicit 4>„ dependence of Recall from Equation B-8 that 


V 


= V 

Ofl,P 


+ V 


C^iyP 


cos<^ + V 

SfL,W 


sin0 

~ o 


(44) 


By applying this and substituting cos0, = [exp(/0,)+exp(-/0,)]/2 and sin<^„ = [exp(/<^,) - 
exp(-/0,)]/2i, the integral becomes 




K L 

OfiyP O 
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(45) 


i 


where 


2x J 0 ' 

h 

2x Jo 


_ f 

2x J 0 


^-ik^S^r^cos{±4>'-4>;) 


d4>„ 


(46) 


With a change of integration variables to </>" = ± </>'-<^„-x/2, these can be converted to standard 
integrals. The result is 


J ^ ^i{mB-\){±4>'-^l2) 




(47) 


Combination of the above equations and substitution of blade coordinates with the radial and 
chordwise counting indices x„'^x^,] lead to the final result for the thickness 

noise integral 
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(48) 


which is the working equation used for far field thickness noise in the computer code. Note that 
the Bessel functions do not depend on the chordwise index p and thus can be computed outside 
this summation. The chordwise sum ZV,^ ,exp[-i(kJ^„'+mB<t>J] corresponds to the chordwise 
Fourier transforms that lead to the functions in earlier formulations. In fact, Equation 48 
can be converted exactly to the earlier forms by setting the inflow angle to zero, placing the 
sources on the advance surface, and changing the summations back to integrals. 


Loading Noise 

The thickness noise integral above has a steady thickness source term plus 2 terms from the 
IP (once-per-revolution) component. The major difference in development of the far field 
loading noise formulas is that all of the P orders appear in an infinite series related to the 
loading harmonics. 


The near field loading noise integral was given by Equation 38 in which the Green’s function 
is the same as in the thickness noise derivation. The derivatives with respect to the source 
spatial variables y, can be written down immediately from Equation 40 


dy^ 


r 30 

o ~ o 

dx'o 


-ikJ,cos(±4>'-4>JG^ 

-ikJ^sm(±4>'-<l>^G„ 


(49) 


In Equation 38, , expresses the time dependence of the loading using tangential angle as the 

variable. To proceed with an analytical form for the tangential integral, the loading must be 
expressed as a Fourier series in 0 q. This is discussed in some detail in Appendix B. If we 
express the loading waveform in terms of loading harmonics according to 

jt=-o& 

as in Equation B-19, then the form required for Equation 38 was Equation B-20 






it = -00 


(51) 
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Substitution into Equation 38 gives 


p ^ pyp 'y 1 j-2'„ 


(52) 


and then application of Equations 49 results in 




-iBk 


47t5„ 






(53) 


where the tangential integrals are 


2tt J 0 ^ 




''‘27jo 


(54) 


With the same change of variables used above in the thickness noise derivation, it can be shown 
that 


j ^ mB -k ^i(„B-kx±4,'-iri2) T fh V 

kSr 


(55) 


where the prime on the Bessel function in the first line denotes the derivative of the function 
with respect to its argument. With these results, the far field loading noise integral becomes 


ikr 


Lm 


-iBk e‘"” 

TTt 

4x5 


52 52 


!i,r k 


Xp5. 


(56) 


where, as with the thickness noise derivation, the blade geometry indices p,v are shown 
explicitly. Again, note that although the Bessel functions are shown under the chordwise 
summation (on p), they are not dependent on this variable. Equation 56 is the working form of 
the loading noise for the far field. 


Like the thickness noise integral. Equation 56 can be changed to earlier forms by setting the 
inflow angle to zero, converting the summations back to integrals, and placing the sources on 
the advance surface. The result is exactly equal to Equation 1 1 in Ref. 9. 

Conclusions on this and other sections are collected in Section 7. 
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SECTION 4 

NUMERICAL PRECISION AND MESH SIZES 


The computer program WOBBLE was exercised to identify the minimum number of radial 
and chordwise mesh points needed to obtain numerical accuracy of 0.5 dR Perfect numencal 
accuracy is defined here as the point where additional mesh points provide no change in the 
calculated noise. The 0.5 dB criterion was chosen because it corresponds to typical 

measurement accuracy. , „ „ 

Parameters that influence the mesh size requirement are blade count, blade geometry, loading 

distribution, tip rotational Mach number, flight Mach number, and harmonic of blade passage 
frequency. The first order effects are blade count, harmonic of BPF, and flight Mach number 
since they influence the wavelength directly. The influence on computational accuracy of other 
parameters is relatively insignificant. Major differences were observed as functions of noise 
harmonic mB and flight speed. Results of the investigation are given in Table 1 on the next 
page for 3 propeller/condition combinations. This table was developed based on single surface 
calculations. If source data are input on both the pressure and suction surfaces of the blades, 

then the tables should be interpreted as applying to each surface. u k 

An 8 bladed, highly swept Prop-Fan was used to generate Table la for a 0.80 Mach number 
cruise condition. From this table, it can be seen that 0.5 dB accuracy can be obtained at blade 
passage frequency, mB = 8, with 10 radial and 5 chordwise grid points. For 2xBPF, mB - 
16 10 radial points and 10 chordwise points are needed for 0.5 dB accuracy. For harmonics 
3 4 and 5 0 5 dB accuracy can be obtained with 20 radial and 15 chordwise points for 3xBPb, 
40 rkdial aid 20 chordwise for 4xBPF, and 40 radial and 25 chordwise points for 5xBPF (mB 

Mmiy fewer grid points are needed for slower commuter propeller applications. Table lb 
gives the numbers of radial and chordwise positions needed to obtain 0.5 dB accuracy for a 
straight four bladed commuter propeller at a flight Mach number of 0.50. It can be seen that 
11 harmonics of BPF (mB = 44) can be obtained with as few as 10 radial and 5 chordwise 
positions while 20 harmonics (mB = 80) can be obtained with 30 radial and 10 chordwise 
positions. Similar results can be seen in Table Ic for the same 4 bladed propeller at 0.20 Mach 

"^Th^ results summarized above are to be used as guidelines only. For any specific propc^^cr 
and/or flight speed an optimum number of radial and chordwise positions must be selected by 
trial and error to obtain the fastest running time. Conservative users can use 50 radial and 40 
chordwise positions. This will provide good accuracy up to 0.80 Mach number and 40 harmonics 
of the shaft rotational frequency although the running time may be unnecessarily long. 
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TABLE 1 

Maximum Harmonic Order {mB) Obtainable at Numerical Accuracy of ±0.5 dB 
As Function of Numbers of Chordwise (Jtv) and Spanwise (#/z) Load Points 

a) Highly Swept Prop-Fan at 0.8 and 800 fps Tip Speed 
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b) Straight Blade Propeller at 0.2 and 800 fps Tip Speed 
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c) Straight Blade Propeller at 0.5 M, and 800 Ips Tip Speed 
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SECTION 5 

PARAMETRIC STUDIES AND THEORETICAL TRENDS 

The radiation equations just derived include several new features that are explored this 
section. The effect on computational precision of accurate source placement is evaluated first. 
Then some source effects are investigated to establish their importance. The third item is the 
main subject of this report, the effect of variation in inflow angle. This is studied in some detail 
using circumferential and axial directivity patterns in both the near and far fields. 


PRECISION OF SOURCE PLACEMENT 

The theory allows the thickness and loading source terms to be placed on the actual blade 
surface (for best precision) or on another reference surface, if desired. The baseline case used 
for comparisons was a full-surface computation for the SR7L Prop-Fan at the design cruise 
condition. The on-surface loads were obtained from an unsteady Euler analysis at the cruise 
design point (Ref 18). Computed axial directivity is shown in Figure 5. Figure 6 shows the 
effect of collapsing loads and thickness source terms onto the camber, chordal, and helicoidal 
advance surfaces. It can be seen that, when the sources were placed on the reference surfaces, 
the errors were up to 1 dB at BPF, 1.5 dB at 2xBPF and 2 dB at 3xBPF. It was expect^ that 
the largest deviation from the on-surface calculation would be the helicoidal surface while the 
least deviation would be the camber line, as the camber line is the closest to the actual surface. 
This does not appear to be the case as each surface can only be classified as being different 
from the on-surface calculation. It appears that the largest difference was due to the collapsing 
of the upper and lower surfaces onto a single zero thickness surface while the differences 
between the helicoidal, chordal, and camber surfaces were secondary. 

For lifting line blade definition and loads, the chordal surface definition is convenient. 
When this is used. Figure 6 indicates that accuracy within 1.5 dB can be expected. 


SOURCE STRENGTH EFFECTS 

Angular inflow produces only one new source effect, the unsteady thickness component. 
This and two other effects, profile drag and radial loading, that were previously known but not 
easily quantified by other methods are also evaluated here. 

Unsteady Thickness Noise 

In Equations 37 and B-8 it was shown that the thickness source can be expressed as the sum 
of a steady term and terms that vary sinusoidally due to the varying blade section velocity. 
These terms were "turned on and off to evaluate their importance for a typical 4 bladed 
commuter propeller flying at 0.50 Mach number. The calculation was made with the propeller 
axis pitched up by 20° to exaggerate the effect. Results are shown in Figure 7 at 90 to the 
flight path since the thickness noise tends to maximize close to the plane-of-rotation. For the 
base case (0° tilt) the source and relative velocity are steady and there is no circumferential 
variation. The solid curve shows the effect of tilting the axis but without the unsteady source 
terms in Equation B-8. The circumferential directivity is due to variations in radiation 
efficiency: for the observer at 90° (under the propeller), the Mach number of the source relative 
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to the observer is higher than for the observer at 270°. The dashed curve shows the effect of 
the unsteady source terms. These add about 2 dB to the peak noise obsen^ed under the 
propeller. 

Radial Loading Noise 

The loading noise integral in Equation 28 or 38 includes three terms via the implied sum 
over the index L The index takes on the values 1,2, and 3, which correspond to the radial, 
tangential, and axial force directions. Tangential and axiaTcomponehts are included in all noise 
theories, but most have neglected the radial component. Although the radial force is generally 
small, there will be a small amount because of blade lean and taper. This was easy to evaluate 
from the Nallasamy load data^'*^ used to generate the baseline case in Figure 5 by simply zeroing 
the radial term. The result shown in Figure 8 indicates that the radial loading component has 
its largest effect 0.25 diameter aft of the propeller plane and less effect forward and aft of that 
point. The maximum contribution of the radial loading component is only about 0.6 dB. 

Note that the loads came from an Euler code which probably does not model the suction due 
to flow accelerating around the tip from the pressure side of the blade to the suction side. This 
effect was discussed in some detail in Ref. 19 where it was shown to be significant under some 
circumstances. 

Effects of Profile Drag 

In lifting line aerodynamic calculations, induced drag and profile drag are treated separately 
in the use of 2D airfoil lift and drag data. For the induced effect, a part of the 2D lift points 
in the direction of section advance due to the 3D induction. Since the profile drag component 
will usually be much smaller than this component, we expect its contribution to noise to be 
small. To verify this expectation, a conventional propeller rather than a Prop-Fan was chosen 
b^ause the blades are relatively thick (5.5 percent of the chord at the 0.75 radius), providing 
higher amounts of profile drag. The power coefficient, advance ratio and Mach number for this 
case are 0.285, 2.17 and 0.50 respectively. It can be seen in Figure 9 that the profile drag 
increases both BPF and 2xBPF by less than 1 dB in the propeller plane and has less effect 
forward and aft of the propeller plane. It is concluded that profile drag has a small to negligible 
effect on the total noise prediction for an on-design condition. (Thickness noise is 10 dB below 
the loading noise levels and doesn’t affect the conclusions.) For very light loading, the relative 
effect could be larger. ’ 


ANGULAR INFLOW EFFECTS 

This section evaluates the magnitude of the angular inflow effects by a series of axial and 
circumferential directivity patterns. 

Axial Directivity and Spectrum 

Calculations in Figure 10 are for a 4 blade propeller operating at a 0.20 Mach number take- 
off condition with advance ratio of 0.895 and power coefficient of 0. 172 . Both the steady and 
once-per-revolution (IP) quasi-steady loads were calculated using a lifting line method. The 
figure shows the noise calculated 4 ways at the peak circumferential location; (1) with no angular 
inflow or unsteady loading, (2) with steady loading only at a 10° angle-of-attack, (3) with 
unsteady loading only but at 0° shaft angle, and finally (4) with all the terms of the new theory. 

It can be seen that the calculations made with unsteady loading increase the noise at BPF 
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and 5xBPF by about 2 dB at the peak noise location. At BPF the increased noise is fairly 
uniform at all locations while at 5xBPF the unsteady loads show increases in noise level only 
aft of the propeller plane. The pure angular inflow effect (no unsteady loads) shows an increase 
of several dB forward of the propeller plane while decreasing aft. At BPF the peak noise is 
increased by about 1 dB but at 5xBPF the peak noise is increased by about 6 dB. It can be seen 
that at BPF both the angular inflow and unsteady loading components contribute significantly 
in modifying the sound pressure level and directivity. However, at 5xBPF the angular inflow 
alone is the major effect. 

Figure 1 1 shows the the effects at 2x, 3x and 4xBPF are intermediate to those at IxBPF and 
5xBPF. 

In general it appears that the angular inflow has a larger effect at the upper harmonics while 
the IP unsteady loading has a larger effect at BPF. This is in agreement with the earlier 
findings of Mani^'\ It should be noted that this generalization may not be true of a more 
complicated flow field where more upper harmonic content exists in the unsteady loading. 
Similar effects were shown in Ref. 3. 

Circumferential Directivity 

A major difference between the angular inflow theory and the previous theories is that the 
propeller axis is not required to be aligned with the flight path. With previous methods, the 
circumferential directivity was always computed for a fixed axial position, since both the plane 
perpendicular to the flight path and the plane perpendicular to the propeller were the same. 
However, with the current method it must be recognized that results can differ significantly in 
the two planes as shown in Figure 12. This is a comparison of the circumferential directivity 
viewed in both planes passing through the propeller center for the SR7L design cruise condition. 
Both curves represent 2D slices through a 3D directivity pattern and it should be recognized that 
neither one gives the whole picture. With this warning, we note that the remainder of the 
circumferential directivity patterns will be viewed in a plane passing through the propeller hub 
and perpendicular to the flight path. 

Shown in Figure 13 is a comparison of the near-field and far-field circumferential 
directivities for a propeller pitched upward. The prediction was made at 0.50 diameter tip 
clearance. Computed levels at the 90° location (below the propeller) were used to normalize 
the two predictions. It can be seen that near-field theory produces an unsymmetrical (left/right) 
pattern that appears to be skewed toward the 0° to 90° quadrant. The far-field theory exhibits 
left/right symmetry and the expected noise maximum directly below the propeller. The near-field 
skewing is due to the observer "feeling" the higher aerodynamic loading on the downgoing side 
of the disc. In the far-field, the observer does not feel this hydrostatic effect but instead only 
senses the radiation effect, which is symmetric and maximizes below the propeller. At smaller 
tip clearances, the effect is more pronounced as shown in Figure 14 for a conventional propeller 
at a 0.50 Mach number cruise condition with a 0.27 diameter tip clearance. The details of the 
condition and method of producing loading are discussed previously in the Profile Drag section. 
This figure also shows the individual effects and combinations of pure angular inflow and 
unsteady loading. It can be seen that maximum noise occurs between 0° and 90° and not 
directly under the propeller as would be expected from far-field theory. 
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SECTION 6 

EXPERIMENTAL VERIFICATION OF ANGULAR INFLOW THEORY 


This section presents a comparison of theory with noise data from a propeller and from 2 
Prop-Fans. One of the Prop-Fans was tested in flight with microphones on a boom attached to 
the wing outboard of the Prop-Fan; the other data were acquired in wind tunnels under 
controlled conditions with angular inflow. 


PROPELLER DATA COMPARISON 

Propeller data were obtained from a joint test by the DFVLR and FAA in the open-jet 
German-Dutch Wind Tunnel described by Dobrzynski et The propeller had 2 blades with 
a diameter of 2.03 meters. Inflow angles ranged from -7.4° to +1.3°, The microphones were 
in the tunnel flow at staggered sideline distances to prevent them from interfering with each 
other aerodynamcically. Since they were in the flow, no corrections were needed for sound 
refraction effects at the jet shear layer. Steady and unsteady aerodynamic loading components 
were calculated with a lifting line method used at Hamilton Standard for aero design and 
structural analysis. 

Figure 15 shows a comparison of measured and predicted noise trends as a function of 
inflow angle and axial directivity at BPF and 5xBPF. The calculations were made at 
measurement locations which would correspond to positions under the flight path. The 
calculations were made at the actual microphone positions rather than at a constant sideline 
distance. It can be seen that the calculated and measured trends agree over the range of observer 
positions and inflow angles. The comparison of experiment and theory can be seen in better 
detail in Figure 16 which shows directivity patterns for 5 harmonics at -1-7.3° inflow angle. 

Figure 17 compares test and theory for a lower power condition with lower blade angle and 
higher tip speed. The comparison includes four harmonics of BPF for a microphone located at 
90° relative to the flight direction. It can be seen that the measurements and predictions agree 
favorably. Figure IB is a comparison of predicted and measured trends with increasing angle 
of attack for a higher tip speed condition. Again, the measurements and predictions agree well. 

The agreement shown in the figures above is very satisfying and indicates that angular 
inflow effects can be predicted well for relatively lightly loaded propellers. 


SR7A PROP-FAN DATA COMPARISON 

SR7A Prop-Fan data were from the Hamilton Standard SR7A model tested in the NASA- 
Lewis 9 by 15 ft anechoic wind tunnel as described by Woodward^*. Inflow angles ranged from 
-10° to -1-15°. A translating microphone was located in the tunnel flow. In the NASA facility, 
when the inflow angle to the drive rig is varied, the distance to the microphone line changes. 
Woodward corrected the data to a constant sideline distance, as reported in Ref. 21 and the 
comparisons are made here on that basis. 

The test case chosen for comparison approximates the takeoff condition. Blade loading for 
takeoff is difficult to predict accurately because the blade sweep and high section lift coefficients 
cause leading edge and tip edge vortices. To evaluate the vortex loading effects, aerodynamic 
loading was computed using 2 independent methods. One of these was an Euler method which 
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was not designed for prediction of tip edge vortices. The other was the Unified Aeroacoustic 
Program (UAAP), reported in Ref. 22. UAAP is a linearized panel code which models vortex 
loading using the suction analogy from delta wing aerodynamics in an attempt to improve noise 
predictions. 

Euler loading for the Prop-Fan was reported by Nallasamy in Ref. 23. Since the actual 
output data set was not available from the aerodynamic calculations, curves were fitted to the 
figures in Nallasamy’s report. The left part of Figure 19 shows the radial distribution of the 
blade section normal force coefficient extracted from the Euler calculations. 

The right part of Figure 19 shows the loading distribution computed with UAAP. The large 
loading peak at the tip is similar to blade surface pressure measurements of highly loaded Prop- 
Fans as reported in Ref. 24. Existence of vortex flow has been verified in flow visualization 
experiments by Vaczy and McCormick in Ref. 25. 

The noise calculations corresponding to the Nallasamy and UAAP loads are presented in 
Figure 20 for a 0° inflow case. Data sets were available for the blade geometry, since this was 
a Hamilton Standard design. Sources were placed on the mean chord line and radial loading was 
not included in the calculation. Since the UAAP loads result in better agreement at the peak, 
they were used for the remainder of the SR7A angular inflow evaluation. 

Figure 21 is a comparison of predicted and measured trends with increasing pitch angle at 
BPF. Measurements and calculations show similar trends. In the forward quadrant, noise 
increases strongly with increasing pitch angle. The predicted curves tend to cross each other in 
the aft quadrant at about 110°. Measurements show similar behavior at 120°. Agreement at 
the far forward and far aft locations is not as good. The wind tunnel environment is particularly 
challenging at these angles for noise measurements and could contribute to the lack of 
agreement. 

In Figure 22 the data and predictions from Figure 21 are plotted in a different form to 
provide a better understanding of trends with inflow angle. The three parts of Figure 22 are 
cross sections of the previous figure at 70°, 90°, and 1 10° with the levels in each part adjusted 
to cross at 0° inflow angle. The good agreement indicates that angular inflow effects can be 
predicted reasonably well with this theory even for a highly loaded Prop-Fan. However, the 
discrepancy that does exist is attributed to inaccurate aerodynamic loading distributions and not 
to the acoustic theory that uses the loading as input. These aerodynamics loads are difficult to 
predict for highly loaded Prop-Fans and more work is needed to model the effects of sweep on 
tip loading. 


SR7L PROP-FAN TEST ASSESSMENT (PTA) FLIGHT DATA COMPARISON 

SR7L Prop-Fan Data were from the Hamilton Standard 9 ft diameter Prop-Fan tested by 
NASA-Lewis at the 0.8 Mach number cruise condition on the PTA aircraft as reported by Little, 
et The Prop-Fan was wing mounted on a nacelle that could be tilted to angles of -3 ° , - 1 ° , 
and -1-2°. The PTA data were selected because they provided high quality measured noise data 
for a range of inflow angles of a large scale Prop-Fan at a high flight Mach number cruise 
condition (Ref. 26). Also, air loads from a three-dimensional unsteady Euler analysis were 
available for use in the noise calculation. The blade loads were obtained in digital form from 
Euler calculations as reported by Nallasamy, Envia, Clark and Groeneweg in Ref. 12. The 
calculations were made on a sideline with a 0.62 Prop-Fan diameter tip clearance for 5 axial 
boom microphone locations ranging from -t-0.50 diameters (forward) of the Prop-Fan to -1.00 
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diameters (aft) of the Prop-Fan. The inflow angles computed were -0.4°, -1-1.6°, and -f-4.6 
degrees which correspond to nacelle tilts of -3° and -1° (down tilt) and -1-2° (up tilt). The 
differences between the inflow angles and the nacelle tilt angle are due to the airplane angle of 
attack and upwash observed at the Prop-Fan plane. The boom diameter ranged from 4. 1 cm at 
the axial location 0.5 Prop-Fan diameters forward of the Prop-Fan to 15.8 cm at 1.0 diameter 
aft of the Prop-Fan. Since these diameters are neither much smaller than nor much larger than 
the wavelengths, scattering at the boom influenced the sound levels. These were corrected to 
free field values using the scattering theory developed by Hanson and Magliozzi*^. No attempt 
was made to account for refraction in the boundary layer of the boom. Correction values given 
in Table 2 were subtracted from measurements to approximate free field values. 


Table 2 


Scattering Corrections for Boom Microphones 



Mic 1 

Mic 2 

Mic 3 

Mic 4 

Mic 5 

BPF 

O.OdB 

l.OdB 

1.4dB 

1.4dB 

0.9dB 

2xBPF 

0.2dB 

3.4dB 

4.1dB 

4.5dB 

5.5dB 

3xBPF 

l.OdB 

4.5dB 

4.6dB 

4.6dB 

4.8dB 


The cruise Mach number, advance ratio and nominal power coefficient were 0.80, 3.12 and 
1.9, respectively. A comparison of measured and predicted data can be seen in Figure 23 for 
the blade passing frequency harmonic as a function of inflow angle and axial location. It can 
be seen that the theory tends to overpredict at the peak axial noise location at the -0.4° inflow 
case, predict the -fl.6° inflow case fairly well and underpredict the peak noise at the +4.6° 
inflow case. The predicted trends close to the Prop-Fan plane show a severe drop in sound 
pressure level with increasing inflow angle. This appears to be due to cancellation occumng 
betweenlhrckhess and loading noise sources. At the -0.4° inflow case, the calculations produce 
a 2 dB reinforcement between thickness and loading noise sources while the +4.6° inflow case 
produces a 2 dB cancellation between thickness and loading noise sources. Also, it can be seen 
that the sound at the most forward microphone position is significantly overpredicted for all three 
inflow conditions. 

Figure 24 are the predicted and measured trends viewed as a function of increasing inflow 
angle for each axial position at BPF. It can be seen that the trends at the 0.0 X/D and -0.25 
X/D locations are predicted to be opposite to those measured. Forward and aft of these 
positions, the predicted noise trends with inflow angle are better but still do not match the 
measured trends accurately. Figures 25 and 26 are similar comparisons at 2xBPF and 3xBPF. 
At 2xBPF, both the predicted and measured trends are in better agreement for all directivities. 
At 3xBPF a reversed trend can be seen at the -0.25 X/D axial location. The good agreement 
in the propeller case and fairly good agreement for Prop-Fan cases in more controlled wind 
tunnel conditions indicates that the flight test environment is producing unexplained loading 
effects which lead to diffences between measurement and prediction. 

To illustrate the strong influence blade loading has on the predicted noise, loading computed 
using a lifting line method was used as input to the calculation. The unsteady loads were 
computed with a quasi-steady lifting line method with flow distortion and angular inflow 
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provided by Lockheed from an airplane panel method. The comparison of predicted and 
measured noise trends with the lifting line method can be seen in Figure 27 for the peak noise 
level at 0.25 X/D axial location. It is clear that the predicted and measured trends are 
significantly improved from the trends shown in Figure 24. Noted that the lifting line method 
results are not shown here to illustrate its superior capability to predict loading over an Eukr 
method but to illustrate the sensitivity the noise prediction method has to loading. It appears that 
more work is needed to obtain the detailed loading information required for accurate Prop-Fan 

noise predictions. 
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SECTION 7 

CONCLUDING REMARKS 


A new frequency domain noise prediction method for propeller noise has been developed 
based on Goldstein’s convected medium version of the acoustic analogy. Because there are no 
geometrical approximations and no limitations on inflow angle, the theory applies equally to 
propellers, Prop-Fans, and helicopter rotors. The computational formulas are simpler than those 
in time domain methods because there is no need to compute retarded time blade locations. 
They are also simpler than those in earlier frequency domain methods because no Bessel 
functions are needed. The radiation formulas are so simple that they only require about 60 lines 
of code for thickness and loading. Furthermore, low harmonic orders can be calculated on a 
personal computer with acceptable run time. 

The theory is exact in the sense that the predictions will be correct if the blade loading and 
geometry are correct. Reflections from other bodies are not included. Also, the theory applies 
only to angular, not distorted, inflow. 

Parametric studies were run with the new method to show behavior of various source 
components. Instructions are given for placing sources on the actual surfaces of the blades (for 
best accuracy) or on reference locations, such as the chordal surface (for simplicity). The 
simpler scheme can cause errors up to 1.5 dB. The unsteady thickness effect associated with 
varying blade speed was investigated. For the case studied, the unsteady thickness source 
contributed only about 2 dB for a 20° inflow angle. For helicopter type operation, the unsteady 
thickness effect would be much larger. Angular inflow influences loading noise via two distinct 
effects. The first is unsteady loading, which has been treated by many theories in the past. 
The second is the unsteady velocity effect, just mentioned with regard to thickness noise. The 
latter has been ignored until recently in propeller noise theories but is shown in this report to 
be more important than unsteady loading (computed with axial inflow). With the new theory, 
there is no reason to neglect either of these effects. 

Experimental verification was provided by comparisons with propeller and Prop-Fan data. 
Agreement with the propeller data was good over a wide range of inflow angle for 5 harmonics. 

Comparisons with Prop-Fan data for both takeoff and cruise conditions were not as good. 
Since the noise theory is exact, it is believed that the problem lies in obtaining good aerodynamic 
estimates of the blade pressure distributions. The source of these problems may be different for 
the 2 conditions. For the takeoff, vortex flow dominates the tip loading and makes accuracy 
dependent on the ability to estimates vortex loading. Methods to deal with this via the suction 
analogy have been developed in previous work at Hamilton Standard^^\ However, empirical 
factors are involved and more experience is needed to obtain a reliable method. 

For the cruise case, it must be recognized that the blade sweep was designed to produce a 
large amount of noise cancellation from different radial stations on the blades. This noise 
reduction method has been proved to be effective; however, it makes noise prediction more 
difficult because small errors in load estimates produce large errors in noise. 
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APPENDIX A 

PROPELLER COORDINATES, THE GREEN'S FUNCTION, AND ITS DERIVATIVES 


The radiation integrals in the main text of this report are based on the moving medium, 
harmonic Green function, which is a function of source and observer coordinates. The observer 
coordinates are simply x, y, z, a Cartesian system fixed to the aircraft. As shown in Figure 1 , 
these are visual coordinates with x, by definition, aligned with the flight direction; y and z ^e 
arbitrary but will generally be taken with y pointing to the starboard side and z pointing 
down. The basic source coordinates are x„, y„, z„, which coincide with x. y, z. However, 
since the propeller shaft can point in any direction relative to the flight direction, the appropriate 
coordinate transformations from Xg, y„, Zg to the propeller coordinates Tg, 4>g, Xg' must be 
defined so that the Green’s function and its derivatives can be computed. This appendix 
specifies the transformation of propeller source coordinates and the method of specifying inflow 
angles. Also, working forms for the Green’s function and derivatives are listed for reference. 

It was decided to write the transformations with complete generality so that any shaft 
orientation could be specified. If desired, this will allow flight profiles to be specified with 
arbitrary histories of yaw, pitch, and roll of the aircraft. The standard Euler angles (see 
Reference 17, for example) will be used to conform with aircraft flight dynamics conventions. 
In this system, aircraft orientation is given by specifying angles of yaw, pitch, and roll in the 
coordinate system that is fixed to the aircraft. Note that the order of rotations is important; 
yawing 90° and then pitching 90° results in a different orientation than pitching first and then 
yawing. 

We now define transformations for yaw, pitch, and roll and then link these to find the 
desired net transformation. Each rotation is about one of the aircraft-fixed coordinates in the 
direction of a right-hand screw. 




/ „/// 

INTO PAPER 



Xg =XgCO^oiy -y" sin tty 
sinofy +y''cos«y 

// 


(A-1) 


jrfcosa, *z"sina, 

// /// 
yo = Jo 

Zo=-Xg^\r\aj,*Zg(X)Sotp 


(A-2) 


3'f =)'ocosa,.-z'sina, 
zf =)'oSino(, + z'cosQ!, 


(A-3) 
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The desired transformation is obtained by substituting Equations A-2 into A-1 then A-3 into that 
result: 


y o ~ ^2\^o ^21^0 ^ ^ 23^0 

Zg - Cj^Xg + + O33Z0 


(A-4) 


where 


C,, = COSO! cos 

Hyp 

c ,2 = COS sin sin - sin cos (A-5) 

Cj 3 = sin sin + cos ot^ sin cos 


and 


Cj, = sina^cosa^ 

Cj 2 = cos tty cos oi^ + sin sin sin (A-6) 

Cj 3 = sin oi^ sin cos - cos sin 


c,, = -sin« 

Jt p 

C32 = cos sin (A-7) 

c„ = cos a COSO! 

33 p r 


Equation A-4 is the fundamental transformation from the coordinates aligned with the propeller 
shaft (the primed coordinates) to coordinates aligned with the flight direction (unprim^). The 
derivation was based on visual coordinates but also applies to retarded coordinates, as will be 
seen in Appendix D. The inverse transformation is also useful. This can be found by inverting 
Equations A-1 through A-3 and then preceding as above. Or it can simply be written down 
from the property of orthogonal transformations that the matrix of the inverse transformation is 
the transpose of the original transformation: 


y'o = C,2^„ + C22yo + 

~ <^ 13-^0 ^ 233^0 ^ 33^0 


(A-8) 


Note that, throughout this appendix, all distances and the Green’s function are non- 
dimensionalized by the propeller tip radius. The coordinate sets above are all right handed 
regardless of the sense of propeller rotation. However, for the propeller cylindrical coordinates, 
we define 4>o lo be always positive in the direction of rotation as was shown in Figure 3. 
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Thus, 


y', = r^cos4>„ 
z' = ±r„sin<^„ 


(A-9) 


where here, as in the rest of this report, the upper sign applies to right hand rotation and the 
lower sign applies to left hand rotation. When ^nations A-9 are substituted into Equation A-4, 
the relation between propeller cylindrical coordinates (for right hand rotation) and the original 
Cartesian source coordinates aligned with the flight direction is 


= c,X + c,2/*„cos0„ ± 
yo = + c^2r,cos<t>„ ± c^,r^sin<t>^ 

= c,X + Cj^r^cos^^ ± c,,r^sin(j)^ 


(A-10) 


Now, to simplify notation, we define source-to-observer distances 

7 = 3 ; - 3 ,^ (A-11) 

Z = z-z„ 


so that the harmonic Green’s function can be written 


G_ = 


4x5 


(A-12) 


where the wavenumber k„=mBM,. The amplitude and phase radii are 

5 = v^X2 +/?^(1^+Z^) 


and 

(T = i- (MX + 5) 


(A- 13) 


(A- 14) 


Derivatives are denoted by subscripts as in the following 

X, = = dX/dr^ 

x,=x, = (i/oax/a0, 

X 3 = x^, = dx/dx'„ 

X^ = dX/dx^ 


(A-15) 
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Derivatives of the distance deltas are easily found from the above definitions: 

= -Cj2COS(^„ T Cj3$in0^ 

X^ = c,2sin0„ T c,3cos0^ 

^. = -1 


(A-16) 

I 


y, = -C22COS<^„ T C23sin<^>„ 
Y^= C22sin0^ T C23COs<^„ 


(A-17) 


I 


Zr = -C32COS(^^ T C33sin0^ 

■^<(> “ ^32 ^ ^33^08^^ 


(A- 18) 


Derivatives of the amplitude radius are 


S^^[XX^^^\YY^^ZZ)]/S 
5, = [XX,3-^2(yy,+zz^]/5 
S^, = [XX^, + ^HYY^+ zzj)ys 
= -x/s 


Derivatives of the phase radius are 

a, = (MX, + 5,)/^^ 
a , = (MX , + X ,)/|8^ 

= + W 

And, finally, the Green’s function derivatives are 

Gr = (iKa^-S/S)G^ 


(A-20) 


(A-21) 


Equations A-5 through A-7 and A- 10 through A-21 are the working forms for computation of 
the Green’s function and its derivatives for use in the near field radiation equations of the main 
text of this report. 
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APPENDIX B 

COMPUTATION OF SOURCE TERMS FOR RADIATION INTEGRALS 

This appendix provides information on computation of the thickness source V^ , for the 
near field and far field radiation formulas (Equation 37 and 48), the loading source , for the 
near field formula Equation 38, and the loading source harmonics F*^ , for either near or far 
field use. The thickness source term is first given without approximation in a form that can be 
used in full surface application of the radiation integrals. In this application, locations of the 
surface area elements in terms of radial coordinate r^, axial coordinate „ and tangential 
coordinate 4> as defined in Figure 4 must be given. Also, the unit normal vectors must be 
available. For more approximate work, a mean surface form of the thickness source is given 
in terms of airfoil thickness distribution in a manner similar to earlier Hamilton Standard 
formulations. The loading source is also treated in full surface form, which requires unit 
normals and pressure distributions, and in various mean surface forms involving section thrust 
and power or lift and drag. 

Thickness Source - Exact Form 

As defined in the main text, the thickness source is 


where = V^/c„ and, as defined in conjunction with Figure 1, V^r is the velocity of the ar^ 
element normal to itself relative to the fluid and AA is the elemental area normalized by r, . 
More specifically, the normal velocity is 

V. - V . • « (®'2) 


where is simply the velocity of the area element relative to the fluid. 

Motion relative to the fluid arises from 2 elements: forward flight and blade rotation. Recall 
from Figure 3 that the speed K = is in the flight direction whereas the rotation is 
about the propeller axis, x^. It is most convenient to perform the dot product of EquaUon B-2 
in propeller coordinates where it is assumed that the area-element unit normals are available in 
cylindrical coordinates. 

n=(n^,n^,nj 


where points in the direction of rotation. To express in propeller coordinates, note that 
Equation A-8 gave the transformation from flight coordinates to propeller coordinates in terms 
of a matrix with elements Cy. Since this transformation also applies to velocity components, 
the components of can be written down immediately in Cartesian coordinates 


V =c V 

V ‘'ll 

y' 12 X 

V =c V 

V ^13 ^x 


(B-4) 


Figure 28 can be used to convert to cylindrical coordinates. When the rotational velocity 
component is added, we find 
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V;= ± 

K ± V^,cos<l)^ +Ur^ 

V,= V, 


(B-5) 


where, as throughout this report, the upper sign applies to right hand ro^tatfon. Substitution of 
Equations B-4 into B-5 gives the relative velocity comp>onents in prop>eller coordinates. 


F = c,2 V cos ± c,3 V sin 

^■t> " ± CijV COS0^ + Qr^ 

V.= c„K 


(B-6) 


This can now be dotted into the unit normal of Equation B-3 to give, after rearrangement, 

(B-7) 




+ ± c,3«^ V cos0^ + (±c,3«^ - V^sin<i>^ 


Now, by defining the tip rotational and flight Mach numbers as M, = Qr /c^ and M = 
VJc„, redefining r„ as radius ratio, and recalling the definition in Equation B-1, the final form 
can be written as 


V 

A*." 


= V + 

OfJL^W 




+ V sin <6 

StA,9 


(B-8) 


where the components are 

K,., = [(c,2«. ± c,jn^)MAAl / (B-9) 

K... = l(±c,,n^-c„n^MAAj;, 


Equations B-8 and B-9 are the working forms for the exact thickness source in terms of the 
flight Mach number, the rotational Mach number, inflow angle (leading to the c^^’s), and the unit 
normals. Equation B-8 includes a constant term (independent of <f>J corresponding to the usual 
steady thickness noise and 2 terms dependent on which produce a once-per-revolution 
unsteady thickness source component. 


Thickness Source - Mean Surface Approximation 

A convenient alternative to the exact source representation of Equation B-9 can be obtained 
in terms of the airfoil chordwise thickness distribution rather than the unit normals. Geometry 
in the derivation is referred to the camber surface, a refinement over the advance surface used 
in previous formulations. (As will be seen, the same derivation can be interpreted as referring 
to the chordal surface.) The radial component of the surface normal is neglected and 
normals for the upper and lower airfoil surfaces are denoted and By combining 
corresponding elements on the upper and lower surfaces. Equation B-9 becomes 
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K,,. =K±^,iK-^^2^)^^^,.^ 


(B- 10 ) 


Now, for the surface integration in the radiation integrals, the chordwise index »/ runs only 
from leading edge to trailing edge rather than wrapping around the airfoil as with the full surface 
(exact) model. 

Approximate forms are needed for the sums of the unit normal components in Equation B- 
10 . These are deduced from the sketch in Figure 29 with the following result. 

= -sin(0^-0,) +sin(0^+^2) (B-11) 

^u^^2x= cos(fl^ - 0,) - cos(0^ + 


Application of standard trigonometric identities leads to 

^h'^^24. ^ ("COS0, +cos02)sin0^ + (sin0, +sin02)cos0^ 
n^ + ii2x= ( sin 0, + sin 62) sin + (cos0, -co^B^cosO^ 


If the reference surface is the chordal surface, 0 , and B^ are small and the difference of cosines 
is negligible. (If the reference surface is the camber surface, 0 , and B2 are equal and the 
difference of cosines is 0 ). In any case, the sum of the sines is approximately equal to the sum 
of the angles. Thus, we set 

Now, we express the area element as fA/4J^ ,={'ArJ^(A7), where A7 is distance along the 
camber line (normalized by r,). Then, noting that (^1+62) “ tan(0,+^2) dh/dy, it follows 
that (Bi+B2)Ay=Ah, and we can write 

= (r,M,cos sin B)^ (Ah\,, (Ar)^ 

= M(±c, 3cos0^)^ (A/i)^, (Ar)^ (B- 14 ) 

K,,, = ^i2COS0,)^ (A/2)^ , (Ar\ 

Here, B^ is the angle to the camber line (from the plane of rotation) at radius ratio and Ah 
is the change in thickness (normalized by r,) of the panel or area element from its leading edge 
to its trailing edge. Equation B -14 with Equation B-8 are the working forms for the thickness 
source using the camber surface representation. For most cases it will be adequate to use the 
chordal surface for reference. In this case. Equation B -14 applies with B^ set equal to the blade 
angle Ob- 

Before leaving this section, it is worthwhile to relate Equation B -14 to source representations 
used in earlier work at Hamilton Standard. To do this, first note that for axial inflow that C|2 
and Ci3 are 0 so that , and , are 0 . Since c,, =1, then 
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^ ^ ^sin0^) Ah Ar 


(B-15) 


In the earlier approximations, 6^ was the advance angle so that the Mach numbers and advance 
angles were simply related through rJ4, = M^cosB, and M = M.sinO,. Substitution of these 
relations into Equation B-15 leads to , = M^AhAr, which is equivalent to the earlier 
thickness source formulation. 


Loading Source - Exact Form 

Interpretation of the loading source for use in Equation 38 is very easy to state: , is 

simply the force on area element fi,v normalized by ambient pressure and tip radius squared. 
Subscripts / =1, 2, 3 correspond to radial, tangential, and axial components of force acting 
on the fluid (as opposed to acting on the blade). The tangential component is always positive 
in the direction of rotation and thus is ordinarily a positive number. Axial components 
producing positive thrust will be negative numbers because the load pushes the air back along 
the negative axis. 

If the blade pressure distribution and unit normals are available from a finite element 
performance calculation, the loading elements can be computed directly from this information. 
Goldstein’s paper (ref. 13), which was the reference for Equation 1 of the main text, defined 
the loading source as where p and p„ are the surface and ambient pressure 

and is a viscous term which will be ignored here. (We have changed the sign of this 
definition because we are using the outward surface normals rather than the inward normals of 
Goldstein’s paper J The force on an area element is fdA. When this in normalized by pj^, 
and the area is normalized by r/, the load element for Equation 38 is 

„ \ip-pf)n,AA 


(B-16) 


The time reference for is the moment the reference blade pitch change axis is 

aligned with the yf axis. If it is desired to use a different reference, then the loading 
waveform must be lagged accordingly. Referring to Figure 30, say that the waveform is 

specified as (i.e. determined from some analysis to be) where fir = fir' -1- and 

0ref is the angle measured in the direction of rotation from the yf axis to the new reference 
point. Then since ,(0 t) = F.^(firO it follows that 

(B-17) 

Actually, an additional time shift occurred in deriving Equation 38 so that the necessary quantity 
is 


This is all that is needed to apply Equation 38 if a direct waveform input is used. 
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A drawback to the direct waveform approach is that many data points will be needed to 
specify the load at each of the area elements and these will each have to be interpolated to use 
Equation B-17. Because of this, it may be more desirable to input the loading to the noise 
subroutine in terms of loading harmonics. Then only a few harmonic coefficients can be used 
to specify a fairly complicated waveform. Furthermore, the interpolation problem is avoided 
by shifting the time lag to a phase lag. If the loading waveform is specified in terms of complex 
Fourier coefficients as 




-ikfh' 


(B-19) 


i = -00 


then the form needed for Equation 38 is 


(B-20) 


k^-oo 


Similarly, if the loading waveform is defined in terms of ordinary cosine and sine coefficients 

= F^,.. * E [^'a.,..cos(aT0 + F«,„.sin(Or')] (B-21) 


k=\ 


the form needed for Equation 38 is 

00 

+ E COS (0^ <f>^^ + , sin (<f>^ (B-22) 


*=1 


This equation used to generate the loading waveform (using <f>^ as the time parameter) in terms 
of the loading harmonics. 


Loading Source - Mean Surface Forms 

The form for the loading element given above in Equation B-16 is convenient to use if the 
blade pressure distribution is available from a numerical calculation or from an experiment. 
However, often blade loading must be estimated using less sophisticated methods. For example, 
a lifting line code may output radial distributions of thrust and power (or their coefficients) or 
section lift and drag coefficients may be available. This section discusses use of these 
parameters. Note that proper use of these guarantees that all of the loading that affects 
performance will be included. However, these performance methods do not generally provide 
radial loading. Also, chordwise distribution of loading must be estimated. In the following, 
formulas are given for the case where chordwise loading is uniform; the computer code includes 
a 2-parameter function generator that can be used to approximate actual chordwise distributions. 

Note that in the following only the force magnitudes are given. The time history discussion 
given above for the exact form in Equations B-16 to B-22 applies equally here. 
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Power and thrust coefficients - The load element relating to tangential force can be written 


2/x,f 


[J(tang force)/JrJ 
KPor" 


(B-23) 


where the radius elements are non-dimensional and N„ the number of panels across the chord 
(counted by v), divides the loading equally among the panels. Tangential force can be obtained 
from the definition of power coefficient (in gradient form) 


dCp _ <^(power)/rfr^ 


(B-24) 


By noting that ^/(power) = Qr„r,t/(tang force) and that tip speed can be expressed as -kND = 
c^^, it can be shown that 


dCp TT^ (/(tang force) 


(B-25) 


Combining the above equations and noting that leads to 




dC^ 

i^rN dr 


(B-26) 


A similar derivation for the thrust component gives 


3ft, »» 


47M;Ar„ dC, 
■K^ N dr 


(B-27) 


where the minus sign accounts for the fact that positive thrust coefficient forces the air in the 
negative x direction” 

Tangential force and thrust gradients (dimensional) - An unsteady lifting line code may output 
blade loads in pounds per inch of span. If we write the tangential force per unit span as 
Ib/in, then F. is just a number and the tangential source in Equation B-23 becomes 


2fi,7 


(F^ Ib/in) r, Ar^ 
NPor" 


(B-28) 
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To cause the units to cancel, we define the ambient pressure to be ;?„=P„lb/in^ and the tip 
radius to be r,=i?,in. Then, the tangential load element is 


2/*.- 


F. Ar 

4> o 

N P R, 

¥ O t 


(B-29) 


A similar development gives the corresponding axial loading element 


3^4, F 


NPA, 


(B-30) 


where the axial force per unit span is F^lb/in. The reason for defining the new quantities for 
radius and ambient pressure here is to permit computer output to be used directly. For example, 
R, is the number of inches in the radius and is therefore a pure number. 


Lift and drag coefficients - Lifting line programs define lift coefficient to be lift force/unit span 
normalized by a dynamic pressure times the chord. The direction of the lift force is normal to 
the relative air direction, including the induction, sketched in Figure 31. Thus, the lift force on 
a span section of length r,Ar„ is The relative Mach number usually 

includes a small induced effect. The drag force is the same with Q replaced by Q,. When 
the tangential and axial components of these are found based on the relative air angle ?>, equal 
loading along the chord is assumed, and it is recalled that pjcJ‘=yp„, the results are 


2/*.- 


yMl ^ (Qsin$ + C^cos$) 


(B-31) 




yMjBr 


Ar 

— - (- C, cos$ + C„ sin$ ) 

^ V Z. D / 

F 


(B-32) 


where = b/(2rj is the chord-to-diameter ratio. These permit lifting line output to be used 
directly since Q, Q, and are printed out as is the sum of the advance and induced 
angles. 

Of the mean surface forms of loading input derived above (power and thrust coefficients, 
tangential force and thrust gradients, and lift and drag coefficients), only the latter are presently 
included in the code. The others are listed as examples for the many kinds of loading input that 
may be available. 
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APPENDIX C 

BLADE GEOMETRY FOR MEAN SURFACE REPRESENTATIONS 

Figure 4 in the main text defines the locations of the blade surface elements <f>^, for 
use in evaluating the Green’s function in the near field radiation integrals of Equations 37 and 
38 and the far field integrals of Equations 48 and 56. For full surface applications of the theory, 
it is expected that these surface locations would be obtained from a finite element aerodynamic 
code. However, simpler means are required for use with lifting line programs. This section 
describes one such scheme that is based on arrays versus radius ratio that are used in Hamilton 
Standard program H-444 (available at NASA-Lewis), namely the chord and blade angle plus 
x,y,z coordinates of the sweep line. The method provides the arrays X(n,y) and PHI(^,j/) 
which are the x„', coordinates of source locations along the chord line. Recall that the index 
H counts radial locations. For mean surface use, the index v counts elements in the chordwise 
direction, starting at the leading edge. 

We start by finding the coordinates of the mid-chord line where it intersects the cylindrical 
source definition surface at = constant as shown in Figure 32. The H-444 convention is that 
the X coordinate, which we will call x^4^, coincides with the pitch change axis; Z444 coincides 
with the axis of rotation in the downstream direction; and completes a left hand coordinate 
system, as shown in the figure. H-444 uses arrays defining the "xyz" or mid-chord line by x^^ 
= x^, etc. With these notations, it is clear that the tangential position of the mid-chord line is 
given by 

(C-1) 


and its axial position by 


X — — 7 
me ^xyz 


(C-2) 


The source element panels can now be located from the mid-chord position based on the 
blade angle 0^ and the chord b at the source radius using the sketch in Figure 33. The 
distances in the "unwrapped" cylindrical surface are 




me 



u-0.5 

N 

w 



(C-3) 


and 


X, = X + — sin 0. 

s me ^ f ^ 


v-0.5 

N 

r 



(C-4) 


where N, is the number of elements in the chordwise direction. By defining chord-to-diameter 
ratio Bp = b/2r, and collecting terms, the working forms are found to be 
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(C-5) 


+ T 


1 - 


2p-1 

N 


Ba COS0J, 


and 





2p-1 


Bo sinOj 


The following notation is closer to the actual array definition and includes the 
explicitly. 





j 2 j '-1 

N 

Bp cos 0^ 

w 

- 


and 


X(pi,p) 





2p-1 

N 

9 


Bo sinfl^ 


j <* 


(C-6) 


subscript 


(C-7) 


(C-8) 
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APPENDIX D 

OBSERVER COORDINATE SYSTEMS 

The position of the observer with respect to the propeller can be represented in a variety of 
coordinate systems. These can be based on the visual (current) propeller location or the retarded 
(emission) location. Furthermore, the x axis can be aligned with the flight direction (flight 
coordinates) or with the propeller axis of rotation (propeller coordinates). And finally, for any 
of these combinations, Cartesian, cylindrical, or spherical coordinates can be used. This 
appendix gives the notation for all of these combinations that appear to be useful and shows how 
to transform from one to the other. Also, it derives identities for certain combinations of 
coordinates and Mach number that arise in the transformation of the far field Green’s function 
integral from the visual to the retarded system. 

Coordinate Notations 

This section simply lists the various notations used for observer position and gives relations 
between the Cartesian, cylindrical, and spherical systems. 


x-axis 

Type 

Visual 

Retarded 

Flight 

Direction 

Cartesian 

cylindrical 

spherical 

X. y, z 

X, R, 4> 

Xr, y, z 
Xr. R. 4> 
r. 6. <t> 

Propeller 

Axis 

Cartesian 

cylindrical 

spherical 

r ' V ^ 7 ^ 
“^V > Jv * 

xj. RJ, 4>J 

x', y', z' 
x'.R',4>' 
r. 6', (f>’ 


The original wave equation was set up in the visual-flight system x, y, z. The propeller-flight 
system is rotated with respect to this by the angles of yaw, pitch, and roll that were defined in 


Appendix A. Within any of the systems (i.e. within any of the boxes in the table 
relations between the various types of coordinates are standard: 

above). 

Cartesian cylindrical 

X =x 

y = Rcos<f) 

X =x 

R = \Jy^+z^ 

(D-1) 


z = Rsin(f> 

tan =z!y 


Cartesian *> spherical 

r = \]x^ +y^+z^ 
cos 6 =x/s/x^ +y^+z^ 

X = rcos 9 
y = rsinScoscj) 

(P-2) 


tan (^ = z/y 

z = rsin^sinf^ 


Cylindrical «• spherical 

r = \jx'^+R} 
cosB =x/\/x^+R^ 

X = rcos 6 
R = rsinO 

(D-3) 


II 

II 
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Visual Coordinates for Input to the Green’s Function 

The Green’s function in Appendix A that is needed for near field calculations requires the 
observer position to be specified in terms of Cartesian coordinates in the visual system jr, y, z 
aligned with the flight direction. In many cases, these might be input directly. However other 
options must be considered. For example, say that cylindrical or spherical coordinates are 
available. Then, Equations D-1 or D-2 can be applied to find x, y, and z. If, however, the 
observer coordinates are available in the system aligned with the propeller axis, then coordinate 
systems must be rotated. Appendix A derived the transformations between the Cartesian system 
aligned with the flight direction and with the propeller axis in terms of a matrix of coefficients 
c,- that can be computed from the angles of yaw, pitch, and roll using Equations A-5, A-6, and 
A-7. In analogy to Equations A-4, the same coefficients can also be applied to the visual 
observer coordinates as follows. 


X = c„x' + c^{yi + CjjZv 

y ~ ^21^^ ^ ^ 22 ^^ **" ^ 23 ^*' 
z = CJ^x!, + Cj^y' + Cjjz' 


(D-4) 


If cylindrical coordinates are available in the visual system. Equations D-1, with the appropriate 
primes and subscripts, can be used to modify Equation D-4 as follows. 

X = c^^x' + c,2^'cos</>i + CjyJtisin4>i 
y = C 2 fXy + Cj 2 l?yCOS<^v + C 2 jK'sin</>' 
z = Cjjxi + Cj2^'cos(f>' + C2jR'sin4>(, 


Angle Relationships in Retarded System 

The transformations in Equations A-4 and A-8 can also be applied to the retarded observer 
coordinates as follows. 


= c„x' + 0 , 2 / + c,3z' 

y = C 21 X + ^22^ + ^23^ 

z = Cj,x' + Cj2y' + C33Z' 


(D-6) 


and 


x' = c„x^ + C2,y +C3,z 

y' = C,2X’^ + ^22y ^ *“32^ 

z' = c,3X^ + C23y + C33Z 


Relationships between the spherical coordinate angles in the propeller and flight retarded systems 
can be derived by applying the right hand column of ^nations D-2 with the appropriate 
subscripts and primes to both sides of Equations D-6 and D-7 and then canceling the r s from 
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(D-8) 


both sides. The result are the following useful sets of equations. 

cos 6 = c,, COS0' + c,2 sin cos 0' + c,j sin sin 
sin 6 cos ^ =0^^ cos O' + sin ^'cos 4>' + sin 6' sin <f>' 
sin 0 sin 0 = Cj, cosd' + sin0'cos<^' + C33 sin 0' sin 0' 

[ 

[ 

COS0' = Cj, cosd + Cjj sin0cos0 + C3, sin^sin*^ 
sin 6' cos <f>' = cos 6 + sin flcos <f> + sin 6 sin <l> (D- 9 ) j 

sin 0' sin = c,3 cos0 + 0^3 sin0cos<^ + C33 smdsin<f> ■ 


For spherical coordinates, these equations can be used to change from the retarded propeller 
system to the retarded flight system and vice versa. Of course, the radiation distance r is the 
same in both systems. 


Conversion from Visual to Retarded Coordinates in Flight System 


This section gives several relationships that are useful in transforming equations from visual 
coordinates to retarded coordinates. They are all based on the equation relating visual and 
retarded x's, which is easily derived with reference to Figure 34 : To find the x position in 
the retarded system, consider a wave emitted from the origin of the x^, y, z system. Say this 
wave takes time to propagate to the observer traveling at the ambient sound speed c„. The 
radiation distance is r=cj„. In the same amount of time the propeller will move a distance 
Mcj„=Mr along the x axis. Thus, the visual and retarded axial distances are related by 
Xr=x+Mr. This can be rewritten 


X =x^~ Mr 

= x^~ M\jx^,+y^+z^ 


(D- 10 ) 


The inverse of this is found by solving for x^-x, squaring both sides, and solving for x^ 

X, = -^[x + Mv/jr^+/3^(y^+z^) ] 

= 1[x^M^/?T^^] (D-11) 

0 ^ V M j 

= J_ (jc + MS\ 

o) 

The last line in Equation D -1 1 uses the amplitude radius 5 „ measured from the propeller center 
as defined in Appendix E: 
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This can be written in terms of the radiation distance and the Doppler factor by substituting 
Equation D-10 for x and collecting terms. The result is 

5^ = r(l-Mcos0) (D-13) 

Another expression involving the amplitude radius needed in the main text can be found by a 
similar approach 

-L(Mx+SJ=r (D-14) 
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APPENDIX E 

FAR FIELD GREEN'S FUNCTION 


For observer positions in the far field, the Green’s function simplifies considerably. In 
particular, it exposes the tangential variable in an exponential so that integration over this 
source coordinate can be performed analytically. This appendix derives the far field form of the 
Green’s function and the derivatives needed for the radiation integrals of the main text. Results 
are given in terms of the variables r, 0', S„, S^, and 5^, which depend on observer position. 
These have especially simple interpretations in terms of retarded coordinates but formulas are 
also presented to compute them from visual coordinates. 

In Appendix A the harmonic, moving medium Green’s function was given by Equations 
A-10 through A-14: 


ik a 

G„ = - — 


(E-1) 


in which S is called the amplitude radius 

S = ^lix-xf^n(y-yfHz-zf] 


(E-2) 


where and a is called the phase radius 

M{x-x)^S 
o ~ 


(E-3) 


In passing to the far field, the usual procedure is followed, namely retaining only observer 
coordinates for the amplitude radius but retaining terms to first order in source distance -j- 
observer distance for the phase radius. For the denominator, then, the amplitude radius is 


= \!x^ + /3^(y^ +z^) 
= \]x^ + 0^R} 


Expansion of the phase radius to first order in (source coordinate) -j- (observer distance) gives 

a = (E-5) 


This can be written more compactly by recalling from Appendix D that {Mx+S^!0^ is the 
radiation distance r and {x+MS^I0^ is the retarded value of the axial coordinate x;. 


o = r - 




(E-6) 
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This can be further manipulated by rotating the source coordinates into alignment with the 
propeller axis using Equations A- 10, changing observer coordinates to the spherical system with 
Equations D-2, and recalling from Equation D-13 that = r(l-Mcos0) 


a = r-(c,,cos0 + C2iSin0 cos<^ +Cj, sin0sin</>) 
- (Cj2cos0 + C22sin0 cos</) + C32 sin0 sin</>) 
T (CjjCOS0 + C23 sin 0 cos + C33 sin ^ sin <^) 


1 -McosO 

^qCOS(^„ 

1 -McosO 
1 -McosO 


(E-7) 


where the upper/lower sign applies to right/left hand rotation. The combinations of terms in 
parentheses are exactly in the form given in Equations D-9, allowing simplification to 


a = r - 


cosO' 

1 -McosO ° 


sin 0' cos 0' 
1 -McosO 


r cos<t>„ 

o ~ o 


T 


sin O' sin 4>' 
1 -McosO 




(E-8) 


In this mixed notation, O' and <f)' are spherical coordinate angles to the observer based on 
retarded coordinates aligned with the propeller axis and 0 is the angle to the observer from the 
flight direction. Recall that observer coordinate is always measured in the direction of a 
right hand screw about the prop axis whereas the source angle <^»o is positive in the direction 
of rotation. By defining new notation as follows 

S (E-9) 

1 - McosO 


and 


^ sin0' 

1 -McosO 


(E-10) 


and using an obvious trigonometric identity, the final form for the phase radius for far field use 
becomes 

a = r - S^xl - r^cos{4>' T (E- 1 1) 

which, as desired, is expressed in terms of the propeller coordinates r„. <f>„, and and is thus 
suitable for integration over the source coordinates. The working form for the far field Green’s 
function is thus 


Q _ ^ -it.S,r.cos(0' T (E-12) 

The form of this result, with (</>' T <f>o) in the exponential and the <l>' dependence of 5„ gives 
rise to the "wobbling mode" behavior discussed extensively in Reference 3. 
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A quantity needed for the far field thickness formula derivation is the derivative of the 
Green’s function with respect to the source coordinate in the flight direction dG„/dXg. To obtain 
this, we write the far field version of G„ as 




G_ = 




(E-13) 


This requires da/dx„ which is most easily obtained from Equation E-3. 


ba 


-J_ [M + 


*^0 


] 


But this derivative only appears in an amplitude role so that in the far field, 





) 


(E-14) 


(E-15) 


Thus, the Green’s function derivative is 


dx^ 

o 


~i k y 

1(M + ±)G 

^2 sj " 


(E-16) 


The actual combination of parameters needed for Equation 39 in the far field thickness noise 
derivation is 


dG M+x/S 

M-^ - ikG = -ik [M f + 1] G„ 


dx„ 


= -ik 


0^ 

1 +Mr/5 

o 


(E-17) 


By using Equations D-13 and D-14, the quantity in parentheses can be shown to be simply the 
reciprocal of the Doppler factor (denoted dj in the main text). Thus, 


dG 

M ^ -ik G = 

dx 


“/ 


(E-18) 


as needed for Equation 39. 
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The far field Green’s function in Equation E-12 is expressed very compactly by use of the 

observer location variables r, and X, defined above via a 

flieht coordinates. Since these are not the best coordinates for actual input to the noise 

the remainder of this appendix is devoted to means of computing these variables from retard^ 

and visual observer coordinates aligned with the flight direction. If other systems are needed, 

they can easily be derived from the information in Appendix D. 


Amplitude radius is defined by Equation E-4 in terms of the visual j:; if input is in 
terms of the retarded value ;c„ then Equation D-10 can be used for the conversion. Radiation 
distance r is defined by the standard distance formula in retarded coordinates as 


r = 



+ 



(E-19) 


If input is in terms of visual coordinates, then x, can be found from Equation D-11. For </>', 
the angle around the propeller axis in the retarded system, relationships can be denved from 

Equations D- 1 and D-7 . 


= z'ly' 


(E-20) 


In Equation E-18, y ' and z' can be computed from coordinates aligned with the flight direction 
using Equations D-7. For the right hand side of those equations, y and z can be input 
directly or obtained from cylindrical coordinates using Equation D-1 as in 

y=/?cos<^> (E-21) 

z = Rsin<^ 

Here R is the sideline distance and <j> is angle measured around the flight axis. For S^, we 
can write 

„ _ cosd' ^ rcose' ^ (E-22) 

"“1-MCOS0 r(l-Mcosfl) 


where x' can be computed from the top line of Equation D-7 and was discussed above. 
Finally, S, can be written 

_ sin0' ^ /-sing^ _ R' _ )jy'^ +z'^ (E-23) 

" "" 1 -Mcosd r(\-Mcos$) 

where y' and z' can be obtained from Equations D-7 with input either directly in terms of 

x^, y, and z or with variations based on Equations D-11 and E-22. 
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APPENDIX F 
LIST OF SYMBOLS’ 


b 

c„ 

df 

f 

h 


— blade chord measured on cylindrical cutting surface (see Figures 4 and C-2). 
= speed of sound in ambient fluid 
= elements of transformation matrix. See Appendix A 
= 1-MCOS0, the Doppler factor 

force per unit area on blade surface, component in coordinate direction 
= blade thickness (local), normalized by r, 


= fT 


K 

p 

Po 

p' 


r 


= mBM„ wave number in harmonic Green’s function 

= pressure in the fluid or on the blade surface. With subscript, disturbance pressure 
= ambient pressure 
= pressure disturbance, p-p^. 

= radiation distance from center of propeller, see Appendix D 
= source radius on propeller 
= propeller tip radius 

= radius to /i** row of source elements on blade, part of blade geometry definition 


t 

V/ 

X 

y 


observer time 

disturbance velocity components in quadrupole source term 
observer (field point) coordinates, jc,y,z 
source coordinates, x„, y„, z„ 

— x position of source element part of blade geometry defintion 


x,x',x„x!^ = observer coordinates. See Appendix D 
Xo,xf„ = source coordinates. See Appendix A 

y>y'>yv = observer coordinates. See Appendix D 
yo>y'o — source coordinates. See Appendix A 

= observer coordinates. See Appendix D 
= source coordinates. See Appendix A 


z z! 


A = blade area 

B = number of blades 

D/Dt = convective derivative, see Equation 6 


Note: After Equation 27, all parameters are non-dimensionalized as follows. Lengths 
and Green’s functions by areas by pressure by p„, and time by Q. 
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ikfi,w 




skfx.p 


JJx) 

J'n(x) 


M 


M, 

P 

R.R' 


S 

So 

Sc 

Ss 


T! 

9 




op,P 


Cfl,P 


’ Sfl.w 


= discrete loading element in coordinate direction i, time dependent, Equations 34,35 
= discrete loading element in coordinate direction /, k!^ complex harmonic 

= discrete loading element in coordinate direction /, k'' cosine comp>onent 

= discrete loading element in coordinate direction /, sine component 

= time dependent Green’s function, see Equation 3 

= harmonic Green’s function, see Equation 22 and Appendices A and E 
= Bessel function of order n and argument x 
= derivative of Bessel function with respect to argument 
= flight Mach number 
= tip rotational Mach number 

= with various subscripts, sound harmonic levels 

= sideline (cylindrical coordinate) distance to observer in flight coordinates, propeller 
coordinates 

= amplitude radius, see Equation 4 

= amplitude radius for far field use, see Appendices D and E 
= defined in Equation E-9 
= defined in Equation E-9 
= time limit for acoustic analogy integrals 
= quadrupole source element, defined in Equation 2 
= source volume for quadrupole integral 

= source normal velocity, see Figure 1 and discussion Appiendix B 
= discrete form of thickness source term, see Equation 37 and Appendix B 
= steady part of , 

= part of , that varies as cos<^<, 

= part of , that varies as sin <f>„ 


Olr 

Of, 

7 

5 

Oc 

e,d' 

K’ 


= angle of pitch of propeller shaft with respect to flight coordinate system 
= angle of roll of propeller shaft with respect to flight coordinate system 
= angle of yaw of propeller shaft with respect to flight coordinate system 
= \-M^ 

= ratio of specific heats, 1.4 for air 
= dirac delta function (impulse function) 

= angle to chord line from plane of rotation, i. e. blade angle 
= angle to observer from flight direction, shaft axis. See Appendix D 
= tangential coordinate to observer in retarded system in flight, propeller system 
= cylindrical coordinate to source element it ii,v, part of blade geometry definition 
= reference angle that serves as r=0 point for loading waveforms 
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p = fluid density 

p„ = ambient density 

p' = density disturbance, p-p^ 

a = phase distance, see Equation 5 

T = source time 

0 = rotational speed of propeller shaft 




Subscripts 

b 

i = 

j 

k 

m = 

n = 

r = 

s = 


index for counting blades 

index for loading and quadrupole components (denotes Cartesian or cylindrical 
components). Note, when i is a repeated index in any expression, a sum over /= 1 , 
2, and 3 is implied. Can denote r,<f>,x' components 

the second index for the quadrupole components. See comment on i regarding 

repeated indices 

index for loading harmonics 

harmonic of blade passing frequency 

harmonic of shaft rotation frequency 

retarded (or emission) position, see Appendix D 

location of source element 


N = component of Mach number or velocity of element normal to itself 
L = loading component of noise 

Q = quadrupole component of noise 

T = thickness component of noise 

p = index for counting source elements in radial direction 

j' = index for counting source elements in chordwise direction 


I 

I 
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APPENDIX G 

USER'S MANUAL FOR COMPUTER PROGRAM 


Input Description for the WOBBLE Computer Program 

The WOBBLE user information is provided in two sections. The first section defines the Input 
Control Data. The second section defines the Blade Geometry and Load Data inputs. 


Input Control Data 

The Input Control Data is read by a location-specified, free field input routine which reads data 
between columns 1 and 72. The first record read in the data set is a description title card that 
can have up to 80 characters. Data begins on the next record. Locations for input are indicated 
by an "L" and followed by a number indicating the desired start of a location. Comment cards 
can be used in this routine to enable quick recognition of input variables. Comments must be 
defined by starting the line with a "C" in column 1. The Input Control Data is terminated by 
an "E" in column 1 . 

Cases can be stacked to run sequentially. After each case the Input Control Data is zeroed and 
must be fully reset for the next case. The run is completed when the word "QUIT" is found on 
a line where a title card is expected. This routine does not allow for data to be input using 
scientific notation (e.g. 0.59E-I-02 is not allowed). Integer input is acceptable. 

An example of a data set can be seen below. Note that blank records can be used between 
cases, if desired. 


(starting in column 1) 

THIS IS THE TITLE FOR CASE 1 
C this is an optional comment 
L 1 1.1 2 3.3 

C this is an optional comment 
Lll 1.0 1 0. -1 

E 

THIS IS THE TITLE FOR CASE 2 
C this is an optional comment 
L 1 7.1 7.2 7.3 

C this is an optional comment 
Lll 0.0 1.0 1.0 1.0 

E 

QUIT 
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This data set will 

cause the WOBBLE program to execute the first case with the 


following values : 



Location 

Value 


1 

1.1 


2 

2.0 


3 

3.3 


11 

1. 


12 

1. 

i 

1 

13 

0. 

t 

14 

-1. 

1 

The second then executes with these values : 

* 

Location 

Value 

: 

1 

7.1 

: 

2 

7.2 

3 

3 

7.3 

i 

11 

0. 

1 

12 

1. 

1 

f 

13 

1. 

1 

14 

1. 

1 

The RUN will then terminate because the word "QUIT" is found where the next title card is 

i 

expected. 


1 

The Input Control Data is always read from UNIT 5. The specific input description is 

1 

as follows : 


1 

input Location 

Description 

t 

1 

Yaw angle (degrees). 


2 

Pitch angle (degrees). 


3 

Roll angle (degrees). 

1 

4 

Input unit number for blade geometry and loads. (Defaults to 

i 

1 


Unit 5) 

- 

5 

Output unit number. (Defaults to Unit 6) 

1 


11 


Type of input indicator: 

1 for full surface, 0 for mean surface. 


12 


Far-field/near- field indicator: 

1 for far-field, 0 for near-field. 


13 


Observer field point input indicator: 

1 for retarded positions, 0 for visual positions. 
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Input Location 

Description 

14 

Propeller direction of rotation indicator: 

1 for right hand rotation, -1 for left hand rotation. (Defaults to right 
hand rotation). 

15 

Unsteady load reference angle (degrees). (Explained in Appendix 
B). 

16 

Type of airfoils: 

1 for Biconvex parabolic, 0 for NACA series 16. (Used in mean 
surface thickness noise calculation). 

17 

Observer sideline input indicator: 

0 for observer coordinates aligned with prop axis, 1 for observer 
coordinates aligned with flight axis. 

20 

Number of observer circumferential positions. (Maximum of 20). 
(Defaults to 1). 

21-40 

Observer circumferential positions (degrees). 

49 

Observer sideline distance in units consistent with those used in the 
propeller definition. 

50 

Number of observer axial positions. (Maximum of 20). (Defaults 
to 1). 

51-70 

Observer axial positions in units consistent with those used in the 
propeller definition. 

80 

Number of input radial positions on the blade. (Maximum of 50). 

81 

Number of input chordwise positions on the blade. (Maximum 
of 80). (Used only if input location 11 = 1.) 

85 

Starting harmonic of Blade Passage Frequency to be computed. 

86 

Ending harmonic of Blade Passage Frequency to be computed. 
(Maximum of 20). 

88 

Number of input load harmonics. (Maximum of 10). 

90 

Number of azimuthal integration steps/ mB. (Defaults to 8). Number 
of azimuthal integration steps are increased as a function of 
harmonic of BPF as follows : 

Number of steps = (Location 90) * mB, where mB is the harmonic 
of Blade Passage Frequency multiplied by the Number of Blades. 
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Input Location Description 

91 Minimum number of radial integration steps. (Defaults to 10). 

(Applies to Mean Surface input only). 

r 

92 Minimum number of chordwise integration steps. (Defaults to 10). | 

(Applies to Mean Surface input only). ' 

93 Increment in number of radial integration steps. Increases the i 

number of radial calculations at each harmonic of Blade Passage I 

Frequency. (Defaults to 0). (Applies to Mean Surface input only). * 

94 Increment in number of chordwise integration steps. Increases the | 

number of chordwise calculations at each harmonic of Blade | 

Passage Frequency. (Defaults to 0). (Applies to Mean Surface : 

input only.) t 

NOTE: Locations 93 and 94 are provided to allow the user to compute at finer grid on 

the blade surface at increasing BPF harmonics. The following conditions must be ^ 

met for array subscripts not to be exceeded: ! 

(Location 91) + (Location 86) * (Location 93) < = 50 j 

(Location 92) + (Location 86) * (Location 94) < = 80 | 

I 

Blade Geometry and Load Data | 

I 

The Blade Geometry and Load Data are read from the unit defined by Location 4 of the 
Input Control Data. This unit can be the same as that of the Input Control Data (Unit 5). If 
the unit number is the same as that of the InputTLonirol Data, the Blade Geometry and Load 
Data are appended to the end of each case found in the Input Control Data. This data should 
be placed after the "E" (case end card) found in the Input Control Data, but before the next 
TITLE DESCRIPTION CARD. 

All inputs defined below for the Mean Surface and Exact forms are read in 
"Free-Format". 

Mu is the number of data points read as a function of radius. The number of radial 
inputs must equal the number shown in Location 80 of the Input Control Data. 

Nu is the number of data points r^d as a function of chord. The number of chordwise 
inputs must equal the number shown in Location 81 of the Input Control Data. 

K is the number of points read as a function of load harmonic. The number of load 
harmonics must equal the number shown in Location 88 of the Input Control Data. All unsteady 
load harmonics are assumed to be based on a single-sided Fourier transform, i.e. positive 
frequencies only. 
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Data are always read in the following order : 

DO 20 k = 1 , Location 88 
DO 10 nu = 1, Location 81 

READ (UNIT,*) (DATA(mu,nu,k), mu = l, Location 80) 

10 continue 
20 continue 

where DATA can be r/R, b/D, Clk, etc. Only the applicable loops (mu, nu, or k) are read. 
The variable dependence on mu, nu, and k is given in the Mean Surface Input and Exact Input 
description sections below. 

Mean Surface Input - This input type is expected when Location 1 1 of the Input Control Data 
= 0. See the first case of sample input for an example of this type. 

The physical description of the Mean Surface Input data is defined in Appendices B and 
C of this report. The WOBBLE program will interpolate the Mean Surface Input data to a mesh 
defined by Locations 91 to 94 of the Control Input Data. 

*** Mx Mt SCO D B P/PO 
*** r/R(mu) 

*** h/b(mu) 

*** b/D(mu) 

*** Mr(mu) 

*** Cl(mu) 

*** Cd(mu) 

*** T444(mu) 

*** Q444(mu) 

*** X444(mu) 

*** Y444(mu) 

Z444(mu) 

*** Gs(mu) 

*** Es(mu) 

The following information is read only if Location 88 of the Input Control Data is > 0. 

*** Gu(mu) 

*** Eu(mu) 

*** Clk(mu,k) 


*** Indicates that the input must be started on a new data record. 
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The variable’s dependence on mu, nu, and k are shown in the description as (mu), (nu), 
(k), (mu,nu,k), etc., and 


Mx - 

Flight Mach Number. 

Mt 

Tip Rotational Mach Number. 

SCO - 

Spinner Radius Normalized by Tip Radius. 

D 

Diameter. 

B 

Number of Blades. 

P/PO - 

Pressure at condition normalized by the pressure at sea level and 59 
degrees F. T 

r/R - 

Input radius ratio (root-to-tip). 

h/b - 

Section thickness-to-chord ratio. 


Section chord-to-diameter ratio. ' 

Section relative Mach Number (with induced effects). 

Section steady lift coefficient. 

Section steady drag coefficient. 

Section operating blade angle (degrees). 

Section air angle (with induced effects, in degrees). 

"X" position of blade mid-chord (see Appendix C). 

"Y" position of blade mid-chord (see Appendix C). 

"Z" position of blade mid-chord (see Appendix C). 

Used in defining the steady chordwise load distribution such that the steady 
chordwise loading shape = Xc * exp(Gs-l) * (1-Xc) * exp(Es-l). Xc is 
the chordwise location, where Xc=0 is the leading edge and Xc=l is 
trailing edge. 

Used in defining the unsteady chordwise load distribution such that the 
unsteady chordwise loading shape = Xc * exp(Gu-l) * (1-Xc) * exp(Eu-l). 

The loading functions are normalized such that the chordwise summation 
of the discrete points equal 1.0. G = 1 and E=1 provide equal chordwise 
loads at each chordwise position. G=2 and E =2 provide a biconvex 
parabolic loading shape while G = 1 and E=5 provide high loading towards 
the leading edge. This function can be used to model chordwise loading 
effects based on blade surface pressure test data, full surface calculations, 
etc. If the actual chordwise loading shape is not known, a constant 
distribution (G = l, E=l) may be used as a reasonable assumption. 

Section unsteady lift coefficient (for k = 1 to 10). Real and imaginary 
parts to be input as a complex number (real, imaginary). Each load 
harmonic must be started on a new record. 

Exact Input Definition - This type of input is expected when Location II of the Input Control 
Data = I . See the second case of sample input for an example of this type. 

Selecting the Exact Input option causes the WOBBLE program to compute the noise for 
only the positions defined by the input. The WOBBLE program will not interpolate this input 
to a finer grid. Tbe physical description of the Exact Input data is defined in the Main section 
and Appendix B of this report. 


b/D - 
Mr 
Cl 
Cd 

T444 - 
Q444 - 
X444 - 
Y444 - 
Z444 - 
Gs,Es - 


Gu,Eu - 


Clk - 
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*** Mx Mt P/PO D B 
*** Rs(mu) 

*** Xs(mu,nu) 

*** PHIs(mu,nu) 

*** Vo(mu,nu) 

*** Vc(mu,nu) 

*** Vs(mu,nu) 

♦** FR(mu,nu,k) 

*** FT(mu,nu,k) 

*** FX(mu,nu,k) 

*** Indicates that the input must be started on a new data record. 


The variable’s dependence on mu, nu, and k are shown in the description as (mu), (nu), 
(k), (mu,nu,k) etc., and 

Mx - Flight Mach Number 

Mt - Tip Rotational Mach Number 

P/PO - Pressure at condition normalized by the Pressure at Sea Level and 

59 Degrees F. 

D - Diameter 

B - Number of Blades 

Rs - Input radius ratio (root-to-tip) 

Xs - "X" location of blade element normalized by (D/2) 

PHIs - Tangential blade element position (in radian) 

Vo - Thickness source term 

Vc - Thickness source term (unsteady term) 

Vs - Thickness source term (unsteady term) 

FR - Radial loading term* (k=0 to 10) 

FT - Tangential loading term* (k=0 to 10) 

FX - Axial loading term* (k=0 to 10) 

* Real and imaginary parts to be input as a complex number (real, imaginary). Each load 
harmonic must be started on a new record. 


Sample Cases for the WOBBLE Computer Program 

Included are sample input and output data sets of three cases that were stacked to be run 
from one execution of WOBBLE. The test cases shown are set-up to be fast running 
verifications that WOBBLE is working correctly on the user’s system. Actual noise calculations 
should be made using a substantially finer mesh to obtain a better blade loading, blade geometry 
and noise radiation representation. 
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Input Data 

The Blade Geometry and Load Data are appended to the Control Input Data so that WOBBLE 
can be run without opening another input unit. 

The first case to be run is a Mean Surface Input case. The Control Input Data is ended 
by "E" and the Blade Geometry and Load Data follow. The first record of data following the 
"E" contains the flight Mach number, tip Mach number, spinner cutoff ratio, number of blades 
and ambient pressure ratio. The next two lines define the input radius ratios. These could have 
been combined on one record, but are split here to demonstrate that the data can spill over to 
the next record. The next record defines the thickness-to-chord ratio, etc. 

The next case in the run is for an Exact Input case. After the Control Input Data ends 
with an "E" the Blade Geometry and Load Data can be seen. The first record following the "E" 
contains flight Mach number, tip Mach number and ambient pressure ratio. The next record 
contains the diameter and number of blades. The next two records contain the radii. The next 
two records contain the Xs corresponding to the radii specified on the previous records for the 
first chordwise position, while the next two records contain the Xs corresponding to the second 
chordwise position for each radius. The remainder of the Exact Input follows the same format. 

The last case in the run is for Mean Surface Input with unsteady loads. The Blade 
Geometry and Load Data are similar to the first case but two unsteady load harmonics have been 
included. It can be seen that the unsteady loads are input in complex format, i.e. 
(real,imaginary). Also the two records preceding the unsteady loads define the chordwise load 
distribution for the unsteady loads. 

The word QUIT following the unsteady loads signals the program to terminate the input. 
Output Data 

The included output data shows the output of the run corresponding to the input data described 
above. It can be seen that the user selected options for each case are given. The yaw, pitch and 
roll angles correspond to input locations 1 to 3, respectively. The type of input (mean -surface) 
was selected by setting input location 11 equal to 1. The number of input radial stations and 
type of airfoils are defined by locations 80 and 16, respectively. The remainder of the user 
selected inputs are similarly tied to the input control data. The propeller input geometry and 
loading are the shown at the 10 input radii. Calculated geometry and force information are then 
given in the form used in WOBBLE’S calculations. The calculated sound pressure levels and 
phases are then given. The phase information is given in the same output format as for the 
sound pressure levels. 

Note that the first and second cases were selected to produce identical results from the two 
different types of input. 
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Sample Input Data 

The following is a sample input data set which will run three noise cases. The input 
format generally follows the guidelines presented above. For these cases, all data is included 
in one file, which will be read from Unit 5. 


WOBBLE Test case (mean surface input *far-field*) 

C INPUT DATA SET FOR WOBBLE 
C 

C YAW PITCH ROLL UNITIN UNITOUT 
L 01 0, 00. 0, 5. 6. 

C 

C lEXACT NEARFR IVORR IDIR PHIREF IBICON 

L 11 0. 1. 0. +1. 0. 0. 

C 

C NAZIMU 
L 20 1. 

C AZIMUTHAL POSITIONS (DEGREES) 

L 21 90. 

C 

C SIDELINE 
L 49 15. 

C 

C NUMX 

L 50 3. 

C OBSERVER LOCATIONS 
L 51 1. 0.0 -1. 

C 

C NUMMUI NUMNUI 

L 80 10. 0. 

C 

C MSTART MEND 
L 85 1. 5. 

C 

C KEND 

L 88 0. 

C 

C MINPHI MINMU MINNU 

L 90 8. 5. 2. 

C 

C INCMU INCNU 
L 93 0. 0. 

E 


.45739 

.67748 

.1540 

13.0 

6. 

.4599 





.1648 

.2110 

.2895 

.3936 

.5140 






.6399 

,7603 

.8644 

.9429 

.9890 






.2704 

.2133 

.1532 

. 1099 

.0748 

.0580 

.0534 

.0492 

.0435 

.0411 

.1008 

.0943 

.0921 

.0921 

.0921 

.0917 

.0846 

.0700 

.0493 

,0333 

.4708 

.4792 

.4976 

. 5292 

.5744 

.6297 

.6882 

.7423 

.7848 

.8098 

-.030 

.030 

. 190 

.384 

.507 

.483 

.526 

,551 

.534 

.417 

.0054 

.0067 

.0068 

.0067 

.0046 

.0041 

.0041 

.0041 

,0044 

.0043 

72.50 

69.77 

65.26 

59.97 

54.09 

48.06 

43.93 

40.86 

38.66 

37.36 

76.12 

72.80 

67.73 

61.59 

55.07 

48.81 

44.12 

40.62 

38.25 

37.75 

.1648 

.2110 

.2895 

.3936 

. 5140 

.6398 

.7599 

.8628 

.9395 

.9839 

.0023 

.0015 • 

-.0001 

-.0022 

-.0017 

.00680 

.0254 

.0520 

,0807 

.1012 

.0170 

.0111 

.0014 

-.0094 

-.0174 

-.0197 

-..0137 

-.0011 

.0128 

.0220 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 

2. 


WOBBLE Test case (exact input version of previous case *far~f ield* ) 
C INPUT DATA SET FOR WOBBLE 
C 

C YAW PITCH ROLL UNITIN UNITOUT 
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00 . 


0 . 


L 

C 

C 

L 

C 

c 

L 

C 

L 

C 

C 

L 

C 

c 

L 

C 

L 

C 

C 

L 

C 

c 

L 

C 

c 

L 

C 

C 

L 

C 

C 

L 

E 


01 0 . 00 . 0 . 

lEXACT NEARER IVORR 

11 1 . 1 . 0 . 


6 . 


IDIR 

1 . 


PHIREF 

0 . 


IBICON 

0 . 


NAZIMU 

20 1 . 

AZIMUTHAL POSITIONS (DEGREES) 
21 90. 

SIDELINE 

49 15. 

NUMX 

50 3. 

OBSERVER LOCATIONS 


51 


1 . 


0.0 


- 1 . 


NUMMUI NUMNUI 
80 5. 2. 

MSTART MEND 
85 1. 5. 

KEND 

88 0 . 

MINPHI MINMU MINNU 


90 8. 0. 

INCMU INCNU 


93 


0 . 


0.457390000000000 
13.0000000000000 
0.238600000000000 
0.746200000000000 
0.357245940077783D- 
0.442964449845996D- 
-0. 51 103702288033 1D■ 
-0.154911293699817D• 
0.689064455673189D- 
0.976495156944241D- 
-0.767643450274055D- 
-0.720153347504015D- 
0.294798633966738D' 
0.106224567292806D- 
-0.288902647357261D 
-0.104100070927512D 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
O.OOOOOOOOOOOQOOO 
0.000000000000000 
0.000000000000000 
( 0 . 000000000000000 , 0 
( 0 . 000000000000000 , 0 
( 0 . 000000000000000,0 
( 0 . 000000000000000,0 
( 0 . 000000000000000,0 
( 0 . 000000000000000,0 
( 0 . 000000000000000,0 
( 0 . 000000000000000,0 
( 0 . 000000000000000 , 0 


0 . 


0. 677480000000000 
6.00000000000000 
0 . 407800000000000 
0.915400000000000 
-001 0.499299471883628D-001 

-001 0 . 100285621154667D-001 

-001 -0.292429039989941D-001 

-001 -0.258897086122820D-001 

-001 0. 629425230930344D-001 

-002 - 0 . 533964461166267D-001 

-001 -0. 524422277106062D-001 

-001 -0.101113516622577 

-002 0 . 176128495258320D-002 

-002 0 . 670000715101688D-003 

-002 - 0 . 172605917030540D-002 

-002 -0.656600669140061D-003 

0 . 000000000000000 
0 . ooooooooooooooo 
0.000000000000000 
0 . ooooooooooooooo 
0 . ooooooooooooooo 
0 . ooooooooooooooo 
0 . ooooooooooooooo 
0 . ooooooooooooooo 

.OOOOOOOOOOOOOOO) 

.ooooooooooooooo) 

.OOOOOOOOOOOOOOO) 

.OOOOOOOOOOOOOOO) 

. OOOOOOOOOOOOOOO) 
.OOOOOOOOOOOOOOO) 
.OOOOOOOOOOOOOOO) 
.OOOOOOOOOOOOOOO) 
.OOOOOOOOOOOOOOO) 


0.459900000000000 
0 . 577000000000000 
0.S42971863335934D-001 
- 0 . 171953594868897D-001 
0.456127159432962D-001 
-0.544639481142936D-001 
0.124183297997186D-002 
-0.121699626169200D-002 
0 . OOOOOOOOOOOOOOO 
0 . OOOOOOOOOOOOOOO 
0.000000000000000 
0 . ooooooooooooooo 
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( 0 . 000000000000000 , 0 . 000000000000000 ) 

(0.2188384196077780-003,0.000000000000000) 

(0.1094236795813350-002,0.000000000000000) 

(0.1548264128289990-002,0.000000000000000) 

(0.1734079936976320-002,0.000000000000000) 

(0.1364695425320520-002,0.000000000000000) 

(0.2188384196077780-003,0.000000000000000) 

(0.1094236795813350-002,0.000000000000000) 

(0.1548264128289990-002 , 0 . 000000000000000 ) 

(0.1734079936976320-002,0.000000000000000) 

(0.1364695425320520-002,0.000000000000000) 

(-0.5666923083533080-004,0.000000000000000) 

(-0.5879961151363220-003,0.000000000000000) 

(-0.1190571994728460-002,0.000000000000000) 

(-0.1726836523758370-002,0.000000000000000) 

(-0.1653403315700340-002,0.000000000000000) 

(-0.5666923083533080-004,0.000000000000000) 

(-0.5879961151363220-003,0.000000000000000) 

(-0.1190571994728460-002,0.000000000000000) 

(-0.1726836523758370-002,0.000000000000000) 

(-0.1653403315700340-002,0.000000000000000) 

WOBBLE Test case (mean surface input with unsteady loads *near-f ield* ) 
C INPUT OATA SET FOR WOBBLE 
C 

C YAW PITCH ROLL UNI TIN UN I TOUT 
L 01 5. 5. 0. 5. 6. 

C 

C lEXACT NEARER IVORR lOIR PHIREF IBICON 

L 11 0. 0. 0. +1. 0. 0. 

C 

C NAZIMU 
L 20 4. 

C AZIMUTHAL POSITIONS (DEGREES) 

L 21 0. 90. 180. 270. 

C 

C SIDELINE 
L 49 15. 

C 

C NUMX 

L 50 1. 

C OBSERVER LOCATIONS 
L 51 0.0 

C 

C NUMMUI NUMNUI 

L 80 10. 0. 

C 

C MSTART MEND 
L 85 1. 1. 

C 

C KEND 

L 88 2. 

C 

C MINPHI MINMU MINNU 

L 90 8. 5. 5. 

C 

C INCMU INCNU 
L 93 5. 5. 

E 


.45739 

.67748 

.1540 

13.0 

6, 

.4599 





.1648 

.2110 

.2895 

.3936 

.5140 

.6399 

.7603 

.8644 

.9429 

.9890 

.2704 

.2133 

.1532 

. 1099 

.0748 

,0580 

.0534 

.0492 

.0435 

.0411 

.1008 

.0943 

.0921 

.0921 

.0921 

.0917 

.0846 

.0700 

.0493 

.0333 

.4708 

.4792 

.4976 

. 5292 

.5744 

.6297 

.6882 

.7423 

.7848 

.8098 

-.030 

.030 

.190 

.384 

. 507 

.483 

. 526 

.551 

.534 

.417 

.0054 

.0067 

.0068 

.0067 

.0046 

.0041 

.0041 

.0041 

.0044 

.0043 
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72*50 

69.77 

65.26 

59.97 54.09 48.06 

43.93 40.86 

38.66 

37,36 

76.12 

72.80 

67.73 

61.59 55.07 48.81 

44.12 40.62 

38.25 

37,75 

.1648 

.2110 

,2895 

.3936 . 

5140 .6398 

.7599 .8628 

.9395 

,9839 

.0023 

.0015 - 

.0001 

-.0022 -. 

0017 .00680 

,0254 .0520 

.0807 

,1012 

.0170 

.0111 

.0014 

-.0094 -. 

0174 -.0197 - 

.0137 -.0011 

.0128 

,0220 

2 , 

2 . 

2 . 

2 . 

2 . 2 . 

2 . 2 . 

2 . 

2 , 

2 . 

2 . 

2 . 

2 . 

2 . 2 . 

2 . 2 . 

2 . 

2 , 

1 . 

1 . 

1 . 

1 . 

1 . 1 . 

1 . 1 . 

1 . 

1 , 

5 . 

5 . 

5 , 

5 . 

5 . 5 . 

5 , 5 . 

5 . 

5 , 

( 0 . 1 , 0 . 

1 ) ( 0.1 

, 0 . 1 ) 

( 0 . 1 , 0 . 1 ) 

( 0 . 1 , 0 . 1 ) (0 

. 1 , 0 . 1 ) 



( 0 . 1 , 0 . 

1 ) ( 0.1 

, 0 . 1 ) 

( 0 . 1 , 0 . 1 ) 

( 0 . 1 , 0 . 1 ) (0 

. 1 , 0 . 1 ) 



( 0 . 01,0 

. 01 ) (0 

. 01 , 0 . 

01 ) ( 0.01 

, 0 . 01 ) ( 0 . 01 , 

0 . 01 ) ( 0 . 01 , 0 . 

01 ) 


( 0 . 01,0 

. 01 ) (0 

. 01 , 0 . 

01 ) ( 0.01 

, 0 . 01 ) ( 0 . 01 , 

0 . 01 ) ( 0 . 01 , 0 . 

01 ) 



QUIT 


Sample Output Data 

The following output resulted from running WOBBLE with the preceding input data. 
The calculation results from the three cases are shown in sequence. It should be expected that 
the general output format will result but differences may appear depending on the output device 
used. The following was printed on a printer that prints 132 characters p>er line and which does 
not recognize the character in column 1 as a printer control character. 


1 

WOBBLE Test case (mean surface input *far-field*) 


SUMMARY OF USER SELECTED INPUTS 


YAW ANGLE = 0.0 DEGREES 

PITCH ANGLE = 0.0 DEGREES 

ROLL ANGLE = 0,0 DEGREES 

TYPE OF INPUT IS MEAN -SURFACE 
NUMBER OF INPUT RADIAL STATIONS = 10 
SERIES 16 AIRFOILS USED IN THICKNESS CALCULATION 


STARTING LOAD HARMONIC <K) = 0 

ENDING LOAD HARMONIC <K) = 0 

REFERENCE PHI = 0.0 DEGREES 


RIGHT HAND PROPELLER ROTATION 


FAR-FIELD THEORY SELECTED 


NUMBER OF INPUT OBSERVER AZIMUTHS = 1 


OBSERVER COORDINATES ARE IN PROP SYSTEM 
(X) INPUTS ARE VISUAL 
NUMBER OF INPUT OBSERVER X(S) = 3 


STARTING NOISE HARMONIC (M) = 1 

ENDING NOISE HARMONIC (M) = 5 
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MINIMUM NUMBER OF TANGENTIAL LOCATIONS USED IN CALCULATION (AT BPF) * 8 

MINIMUM NUMBER OF RADIAL LOCATIONS USED IN CALCULATION (AT BPF) = 5 

INCREASE RADIAL LOCATIONS BY 0*M AT HIGHER HARMONICS OF BPF 
MINIMUM NUMBER OF CHORDWISE LOCATIONS USED IN CALCULATION (AT BPF) = 2 

INCREASE CHORDWISE LOCATIONS BY 0*M AT HIGHER HARMONICS OF BPF 


1 PROPELLER INPUT INFORMATION 


DIAMETER = 13.00 UNIT LENGTH 
NUMBER OF BLADES = 6. 

SPINNER CUT-OFF = .154 


ROR -> 

0.1648 

0.2110 

0.2895 

0.3936 

0.5140 

0.6399 

0.7603 

0.8644 

0.9429 

0 

HB 

0.2704 

0.2133 

0.1532 

0.1099 

0.0748 

0.0580 

0.0534 

0.0492 

0.0435 

0 

BD -> 

0.1008 

0.0943 

0.0921 

0.0921 

0.0921 

0.0917 

0.0846 

0.0700 

0.0493 

0 

MR -> 

0.4708 

0.4792 

0.4976 

0.5292 

0.5744 

0.6297 

0.6882 

0.7423 

0.7848 

0 

CL -> 

-0,0300 

0.0300 

0.1900 

0.3840 

0.5070 

0.4830 

0.5260 

0.5510 

0.5340 

0 

CO -> 

0.0054 

0.0067 

0.0068 

0.0067 

0.0046 

0,0041 

0.0041 

0.0041 

0.0044 

0 

THETA -> 

72.5000 

69.7700 

65.2600 

59.9700 

54.0900 

48.0600 

43.9300 

40.8600 

38.6600 

37 

QPHI -> 

76.1200 

72.8000 

67.7300 

61.5900 

55.0700 

48.8100 

44.1200 

40.6200 

38.2500 

37 

X CHORD -> 

0.1648 

0.2110 

0.2895 

0.3936 

0.5140 

0.6398 

0.7599 

0.8628 

0.9395 

0 

Y CHORD -> 

0.0023 

0.0015 

-0.0001 

-0.0022 

-0.0017 

0.0068 

0.0254 

0.0520 

0.0807 

0 

Z CHORD -> 

0.0170 

0.0111 

0.0014 

-0.0094 

-0.0174 

-0.0197 

-0.0137 

-0.0011 

0.0128 

0 

GAMMA STDY -> 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2 

ETA STDY -> 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2 


FLIGHT, GEOMETERY AND FORCE INFORMATION 


FLIGHT MACH NUMBER = 0.457 
TIP ROTATIONAL MACH NUMBER = 0.677 
PRESSURE / STD SEA LEVEL = 0.460 
PROPELLER DIAMETER (UNIT LENGTH) = 13.00 
NUMBER OF BLADES = 6.00 

CPMEAN = 0.570 

CTMEAN = 0.240 









RO -> 

0.239 

0.408 

0.577 

0.746 0.915 




NU 

z 

1 


XS -> 

0.357E-01 

0.499E-01 

0.543E-01 

0.443E-01 0.100E-01 




NU 

— 

2 


xs -> 

-0.511E-01 

-0.292E-01 

-0.172E-01 

-0.155E-01-0.259E-01 




NU 

_ 

1 


PHIS -> 

0.689E-01 

0.629E-01 

0.456E-01 

0.976E-02-0.534E-01 




NU 

= 

2 


PHIS -> 

-0.768E-01 

-0.524E-01 

-0.545E-01 

-0.720E-01-0.101 




NU 

Jg 

1 


VO -> 

0.295E-02 

0.176E-02 

0.124E-02 

0.106E-02 0.670E-03 




NU 

= 

2 


VO -> 

-0.289E-02 

-0.173E-02 

-0.122E-02 

-0.104E-02-0.657E-03 




NU 

_ 

1 


VC -> 

0.000 

0.000 

0.000 

0.000 0.000 




NU 


2 


VC -> 

0.000 

0.000 

0.000 

0.000 0.000 




NU 

_ 

1 


VS -> 

0.000 

0.000 

0.000 

0.000 0.000 




NU 

z 

2 


vs -> 

0.000 

0.000 

0.000 

0.000 0.000 

K 

_ 

0 

NU 


1 

FR 

REAL -> 

0.000 

0.000 

0.000 

0.000 0.000 

K 

= 

0 

NU 

= 

1 

FR 

IMAG -> 

0.000 

0.000 

0.000 

0.000 0.000 

K 

= 

0 

NU 

= 

2 

FR 

REAL -> 

0.000 

0.000 

0.000 

0.000 0.000 

K 

= 

0 

NU 

= 

2 

FR 

IMAG -> 

0.000 

0.000 

0.000 

0.000 0.000 

K 

_ 

0 

NU 

_ 

1 

FT 

REAL -> 

0.219E-03 

0.109E-02 

0.155E-02 

0.173E-02 0.136E-02 

K 

= 

0 

NU 

= 

1 

FT 

IMAG -> 

0.000 

0.000 

0.000 

0.000 0.000 

K 

s 

0 

NU 

= 

2 

FT 

REAL -> 

0.219E-03 

0.109E-02 

0.155E-02 

0.173E-02 0.136E-02 

K 

= 

0 

NU 

= 

2 

FT 

IMAG -> 

0.000 

0.000 

0.000 

0.000 0.000 


,9890 

,0411 

,0333 

,8098 

4170 

0043 

3600 

7500 

9839 

1012 

0220 

0000 

0000 
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K = 

0 

NU = 

1 

FX 

REAL -> 

-0.567E-04-0.588E-03-0.119E-02-0.173E-02-0.165E-02 

K = 

0 

NU = 

1 

FX 

IMAG -> 

0.000 0.000 0.000 0.000 0.000 

K = 

0 

NU = 

2 

FX 

REAL -> 

•0.567E-04-0.588E-03‘0.119E-02-0.173£-02-0.165E-02 

K = 

0 

NU = 

2 

FX 

IMAG -> 

0.000 0.000 0.000 0.000 0.000 


NU = 1 CHORD(NU)-> 0.250 

NU = 2 CHORD(NU)-> 0.750 


1 SOUND PRESSURE LEVELS IN DB, (RE 20 MICRO PASCALS) 
AT OBSERVER AZIMUTH = 90.0 DEGREES 

AND 15.0 UNIT LENGTH SIDELINE 
PROPELLER COORDINATE SYSTEM 



VISUAL (X) 

1.00 

0.00 

-1.00 


RETARDED (X) 

9.00 

7.72 

6.47 

1X 

BPF 

TOTAL 

118,3 

119.2 

119.8 



MONOPOLE 

106.2 

105.8 

105.0 



DIPOLE 

118,0 

119.0 

119.6 

2X 

BPF 

TOTAL 

107.5 

108.7 

109.2 



MONOPOLE 

100.5 

100.2 

99.4 



DIPOLE 

106.7 

108.1 

108.9 

3X 

BPF 

TOTAL 

98.4 

99.5 

99.9 



MONOPOLE 

94.1 

93.9 

93.0 



DIPOLE 

96.6 

98.2 

99.0 

4X 

BPF 

TOTAL 

89.8 

90.8 

90.9 



MONOPOLE 

87.3 

87.2 

86.2 



DIPOLE 

86.5 

88.4 

89.2 

5X 

BPF 

TOTAL 

81.3 

82.0 

81.8 



MONOPOLE 

79.9 

79.9 

78.8 



DIPOLE 

75.8 

78.0 

78.8 


1 RELATIVE PHASE (IN DEGREES) 

AT OBSERVER AZIMUTH = 90.0 DEGREES 

AND 15.0 UNIT LENGTH SIDELINE 


PROPELLER COORDINATE SYSTEM 



VISUAL (X) 

1.00 

0.00 

-1.00 


RETARDED (X) 

9.00 

7.72 

6.47 

IX 

BPF 

TOTAL 

178.9 

159.1 

142.3 



MONOPOLE 

103.6 

81.3 

62.5 



DIPOLE 

-166.7 

171.4 

152.8 

2X 

BPF 

TOTAL 

104.5 

64.1 

29.7 



MONOPOLE 

39.3 

-5.6 

-43.7 



DIPOLE 

131.0 

86.3 

48.4 

3X 

BPF 

TOTAL 

32.8 

-29.1 

-81.7 



MONOPOLE 

-21.4 

-89.1 

-146.6 



DIPOLE 

70.3 

2.6 

-54.9 

4X 

BPF 

TOTAL 

-40.5 

-124.4 

164.7 



MONOPOLE 

-83.5 

-173.9 

109.3 



DIPOLE 

7.6 

-82.9 

-159.7 

5X 

BPF 

TOTAL 

-115.2 

138.3 

48.7 



MONOPOLE 

-147.4 

99.6 

3.6 



DIPOLE 

-56.9 

-169.9 

94.0 



1 WOBBLE Test case (exact input version of previous case *far-field*) 


SUMMARY OF USER SELECTED INPUTS 


YAW ANGLE = 0.0 DEGREES 

PITCH ANGLE = 0.0 DEGREES 

ROLL ANGLE =: 0.0 DEGREES 

TYPE OF INPUT IS EXACT 

NUMBER OF INPUT RADIAL STATIONS = 5 

NUMBER OF INPUT CHORDWISE STATIONS = 2 


STARTING LOAD HARMONIC (K) = 0 

ENDING LOAD HARMONIC (K) = 0 

REFERENCE PHI = 0.0 DEGREES 


RIGHT HAND PROPELLER ROTATION 


FAR-FIELD THEORY SELECTED 


NUMBER OF INPUT OBSERVER AZIMUTHS = 1 


OBSERVER COORDINATES ARE IN PROP SYSTEM 
(X) INPUTS ARE VISUAL 
NUMBER OF INPUT OBSERVER X(S) = 3 


STARTING NOISE HARMONIC <M) = 1 

ENDING NOISE HARMONIC <M) = 5 

MINIMUM NUMBER OF TANGENTIAL LOCATIONS USED IN CALCULATION (AT BPF) = 8 


1 FLIGHT, GEOMETERY AND FORCE INFORMATION 


FLIGHT MACH NUMBER = 


TIP ROTATIONAL MACH NUMBER = 
PRESSURE / STD SEA LEVEL = 
PROPELLER DIAMETER (UNIT LENGTH) = 
NUMBER OF BLADES = 
CPMEAN = 


CTMEAN = 


0.457 

0.677 

0.460 

13.00 

6.00 

0.570 

0.240 


RO -> 0.239 


0.408 0.577 0.746 0.915 


XS -> 0.357E-01 0.499E-01 0.543E-01 0.443E-01 0.100E-01 

XS ’> -0.511E-01-0.292E-01-0.172E-01-0.155E-01-0.259E-01 


PHIS -> 0.689E-01 0.629E‘01 0.456E-01 0.976E-02-0.534E-01 

PHIS -> -0.768E'01 “0.524E-01-0.545E-01 -0.720E-01 -0. 101 


VO -> 0.295E-02 0.176E-02 0.124E'02 0.106E-02 0.670E-03 

VO -> -0.289E-02-0.173E-02-0.122E-02-0.104E-02-0.657E-03 


NU = 1 

NU = 2 


VC -> 0.000 0.000 

VC -> 0.000 0.000 


0.000 0.000 0.000 

0.000 0.000 0.000 
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NU = 1 

VS -> 

0.000 

0.000 

0.000 

0.000 

0.000 

NU = 2 

vs -> 

0.000 

0.000 

0.000 

0.000 

0.000 


K = 

0 

NU = 

1 

FR 

REAL -> 

0.000 

0.000 

0.000 

0.000 

0.000 

K = 

0 

NU = 

1 

FR 

IMAG -> 

0.000 

0.000 

0.000 

0.000 

0.000 

K = 

0 

NU = 

2 

FR 

REAL -> 

0.000 

0.000 

0.000 

0,000 

0.000 

K = 

0 

NU = 

2 

FR 

IMAG -> 

0.000 

0.000 

0.000 

0,000 

0.000 


1C = 

0 

NU = 

1 

FT 

REAL -> 

0.219E-03 

0.109E-02 

0.155E-02 

0.173E‘02 

0.136E-02 

1C = 

0 

NU = 

1 

FT 

IMAG -> 

0.000 

0.000 

0.000 

0.000 

0,000 

K = 

0 

NU = 

2 

FT 

REAL -> 

0.219E-03 

0.109E‘02 

0.155E-02 

0.173E-02 

0.136E-02 

K - 

0 

NU = 

2 

FT 

IMAG -> 

0.000 

0.000 

0.000 

0.000 

0.000 


K = 

0 

NU = 

1 

FX 

REAL -> -0.567E-04-0.588E-03-0.119E-02-0.173E-02-0.165E-02 

K = 

0 

NU = 

1 

FX 

IMAG -> 0.000 0.000 0.000 0.000 0.000 

K = 

0 

NU = 

2 

FX 

REAL -> -0.567E-04-0.588E-03-0.119E-02-0.173E-02-0.165E-02 

K = 

0 

NU = 

2 

FX 

IMAG -> 0.000 0.000 O.OOO 0.000 0.000 


1 SOUND PRESSURE LEVELS IN DB, (RE 20 MICRO PASCALS) 
AT OBSERVER AZIMUTH = 90.0 DEGREES 


AND 

15.0 UNIT 

LENGTH 

SIDELINE 


PROPELLER COORDINATE SYSTEM 


VISUAL (X) 

1.00 

0.00 

-1.00 

RETARDED (X) 

9.00 

7.72 

6.47 

IX BPF 

TOTAL 

118.3 

119.2 

119.8 


MONOPOLE 

106.2 

105.8 

105.0 


DIPOLE 

118.0 

119.0 

119.6 

2X BPF 

TOTAL 

107.5 

108.7 

109.2 


MONOPOLE 

100.5 

100.2 

99.4 


DIPOLE 

106.7 

108.1 

108.9 

3X BPF 

TOTAL 

98.4 

99.5 

99.9 


MONOPOLE 

94.1 

93.9 

93.0 


DIPOLE 

96.6 

98.2 

99.0 

4X BPF 

TOTAL 

89.8 

90.8 

90.9 


MONOPOLE 

87.3 

87.2 

86.2 


DIPOLE 

86.5 

88.4 

89.2 

5X BPF 

TOTAL 

81.3 

82.0 

81.8 


MONOPOLE 

79.9 

79.9 

78.8 


DIPOLE 

75.8 

78.0 

78.8 


1 RELATIVE PHASE (IN DEGREES) 

AT OBSERVER AZIMUTH = 90.0 DEGREES 

AMD 15.0 UNIT LENGTH SIDELINE 


PROPELLER COORDINATE SYSTEM 


VISUAL (X) 

1.00 

0.00 

-1.00 

RETARDED (X) 

9.00 

7.72 

6.47 

BPF 

TOTAL 

178.9 

159.1 

142.3 


MONOPOLE 

103.6 

81.3 

62.5 


DIPOLE 

-166.7 

171.4 

152.8 

BPF 

TOTAL 

104.5 

64.1 

29.7 


MONOPOLE 

39.3 

-5.6 

-43.7 


DIPOLE 

131.0 

86.3 

48.4 

BPF 

TOTAL 

32.8 

-29.1 

-81.7 


MONOPOLE 

-21.4 

-89.1 

-146.6 


DIPOLE 

70.3 

2.6 

-54.9 

BPF 

TOTAL 

-40.5 

-124.4 

164.7 


MONOPOLE 

-83.5 

-173.9 

109.3 


DIPOLE 

7.6 

-82.9 

-159.7 
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5X BPF 


TOTAL -115,2 138.3 48.7 

MONOPOLE -147.4 99.6 3.6 

DIPOLE -56.9 -169.9 94.0 


1 WOBBLE Test case (mean surface input with unsteady loads *near- field*) 


SUMMARY OF USER SELECTED INPUTS 


YAW ANGLE = 5.0 DEGREES 

PITCH ANGLE = 5.0 DEGREES 

ROLL ANGLE = 0.0 DEGREES 

TYPE OF INPUT IS MEAN-SURFACE 
NUMBER OF INPUT RADIAL STATIONS = 10 
SERIES 16 AIRFOILS USED IN THICKNESS CALCULATION 


STARTING LOAD HARMONIC (K) = 0 

ENDING LOAD HARMONIC (K) = 2 

REFERENCE PHI = 0.0 DEGREES 


RIGHT HAND PROPELLER ROTATION 


NEAR-FIELD THEORY SELECTED 


NUMBER OF INPUT OBSERVER AZIMUTHS = 4 


OBSERVER COORDINATES ARE IN PROP SYSTEM 
(X) INPUTS ARE VISUAL 
NUMBER OF INPUT OBSERVER X(S) = 1 


STARTING NOISE HARMONIC (M) = 1 

ENDING NOISE HARMONIC (M) = 1 


MINIMUM NUMBER OF TANGENTIAL LOCATIONS USED IN CALCULATION (AT BPF) = 8 

MINIMUM NUMBER OF RADIAL LOCATIONS USED IN CALCULATION (AT BPF) = 5 

INCREASE RADIAL LOCATIONS BY 5*M AT HIGHER HARMONICS OF BPF 
MINIMUM NUMBER OF CHORDWISE LOCATIONS USED IN CALCULATION (AT BPF) = 5 

INCREASE CHORDWISE LOCATIONS BY 5*M AT HIGHER HARMONICS OF BPF 


1 PROPELLER INPUT INFORMATION 


DIAMETER = 13.00 UNIT LENGTH 
NUMBER OF BLADES = 6. 

SPINNER CUT-OFF = .154 



ROR -> 

0.1648 

0.2110 

0.2895 

0.3936 

0.5140 

0.6399 

0.7603 

0.8644 

0.9429 

0.9890 


HB -> 

0.2704 

0.2133 

0,1532 

0.1099 

0.0748 

0.0580 

0.0534 

0.0492 

0.0435 

0.0411 


BD -> 

0.1008 

0.0943 

0.0921 

0.0921 

0.0921 

0.0917 

0.0846 

0.0700 

0.0493 

0.0333 


MR -> 

0.4708 

0.4792 

0.4976 

0.5292 

0.5744 

0.6297 

0.6882 

0.7423 

0.7848 

0.8098 


CL -> 

-0.0300 

0.0300 

0.1900 

0.3840 

0.5070 

0.4830 

0.5260 

0.5510 

0.5340 

0.4170 


CD -> 

0.0054 

0.0067 

0.0068 

0.0067 

0.0046 

0.0041 

0.0041 

0.0041 

0.0044 

0.0043 


THETA -> 

72.5000 

69.7700 

65.2600 

59.9700 

54.0900 

48.0600 

43.9300 

40.8600 

38.6600 

37.3600 


QPHI -> 

76.1200 

72.8000 

67.7300 

61.5900 

55.0700 

48.8100 

44.1200 

40.6200 

38.2500 

37.7500 

X 

CHORD -> 

0.1648 

0.2110 

0.2895 

0.3936 

0.5140 

0.6398 

0.7599 

0.8628 

0.9395 

0.9839 

Y 

CHORD -> 

0.0023 

0.0015 

-0.0001 

-0.0022 

-0.0017 

0.0068 

0.0254 

0.0520 

0.0807 

0.1012 

Z 

CHORD -> 

0.0170 

0.0111 

0.0014 

-0.0094 

-0.0174 

-0.0197 

-0.0137 

-0.0011 

0.0128 

0.0220 

GAMMA STDY -> 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

ETA STDY -> 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0D00 

GAMMA 

UNSTD -> 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 
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ETA UNSTD -> 5.0000 5.0000 5.0000 

5.0000 

5.0000 5.0000 5.0000 

5.0000 

5.0000 

1 

CL REAL UNSTD -> 0.100 

0.100 

0.100 

0.100 

0.100 0.100 0.100 

0.100 

0.100 

1 

CL IMAG UNSTD -> 0.100 

0.100 

0.100 

0.100 

0.100 0.100 0.100 

0.100 

0.100 

2 

CL REAL UNSTD -> 0.100E-1 

0.100E-1 

0.100E-1 

0.100E-1 0 

.100E-1 0.100E-1 0.100E-1 0, 

.100E-1 

0.100E-1 

2 

CL IMAG UNSTD -> O.lOOE-1 

0.100E-1 

0.100E-1 

0.100E-1 0 

.100E-1 0.100E-1 0.100E-1 0. 

.100E-1 

0.100E-1 

1 

CD REAL UNSTD -> 0.000 

0.000 

0.000 

0.000 

0.000 0.000 0.000 

0.000 

0.000 

1 

CD IMAG UNSTD -> 0.000 

0.000 

0.000 

0.000 

0.000 0.000 0.000 

0.000 

0.000 

2 

CO REAL UNSTD -> 0.000 

0.000 

0.000 

0.000 

0.000 0.000 0.000 

0.000 

0.000 

2 

CD IMAG UNSTD -> 0.000 

0.000 

0.000 

0.000 

0.000 0.000 0.000 

0.000 

0.000 


1 FLIGHT, GEOMETERY AND FORCE INFORMATION 


FLIGHT MACH NUMBER 
TIP ROTATIONAL MACH NUMBER 
PRESSURE / STD SEA LEVEL 
PROPELLER DIAMETER (UNIT LENGTH) 
NUMBER OF BLADES 
CPMEAN 
CTMEAN 


0.457 

0.677 

0.460 

13.00 

6.00 

0.570 

0.240 


RO -> 0.239 


0.408 0.577 0.746 0.915 


NU = 1 

NU = 2 

NU = 3 

NU = 4 

NU = 5 


XS -> 0.618E-01 0.737E-01 0.757E-01 0.622E-01 0.208E-01 
XS -> 0.270E-01 0.420E-01 0.471E-01 0.383E-01 0.644E-02 
XS -> -0.769E-02 0.103E-01 0.186E-01 0. 144E-01 -0.793E-02 
XS -> -0.424E-01-0.213E-01-0.100E-01-0.951E-02-0.223E-01 
XS -> -0.772E-01-0.530E-01-0.386E-01-0.334E-01-0.367E-01 


NU = 1 

NU - 2 

NU - 3 

NU = 4 

NU = 5 


PHIS -> 0.113 0.976E-01 0.756E-01 0.343E-01-0.391E-01 
PHIS -> 0.543E-01 0.514E-01 0.356E-01 0.159E-02-0.582E‘01 
PHIS -> -0.393E-02 0.525E-02-0.443E-02-0.311E-01-0.773E-01 
PHIS -> -0.622E-01-0.409E-01'0.445E-01-0.638E-01-0.963E-01 
PHIS -> -0.120 -0.871E-01-0.845E-01-0.965E-01-0.115 


NU = 1 

NU = 2 

NU = 3 
NU = 4 
NU = 5 


VO -> 0.228E-02 0.136E-02 0.961E-03 0.823E-03 0.519E-03 
VO -> 0.581E-03 0.348E-03 0.245E-03 0.210E-03 0.133E-03 
VO -> -0.947E-05-0.567E-05-0.400E-05-0.342E-05-0.216E-05 
VO -> -0.799E-03-0.478E-03-0.337E-03-0.289E-03-0.182E-03 
VO -> -0.199E-02-0.119E-02-0.840E-03-0.719E-03-0.454E-03 





NU 


1 


VC -> 

0.698E-04 0.520E-04 0.400E-04 0.344E-04 6.208E-04 




NU 


2 


VC -> 

0.178E-04 0.133E-04 0.102E-04 0.878E-05 0.536e-05 




NU 

- 

3 


VC -> 

-0.290E-06-0.216E-06-0.167E-06-0.143E-06-0.864E-07 




NU 


4 


VC -> 

-0.245E-04-0.182E-04-6.141E-64-0.121E-04-0.728e-05 

i 



NU 

= 

5 


VC -> 

-0.610E -04-0. 455E -04-0. 350E-04-0.301E-04-0.182E-04 




NU 


1 


VS -> 

0.700E-04 0.522E-04 0.402E-04 0.345E-04 0.208E-04 




NU 

- 

2 


VS -> 

0.179E-04 0.133E-04 0.103E-04 0.881E-05 0.532E-05 




NU 

- 

3 


VS -> 

-0.291E-06-0.217E-06-0.167E-06-0.144E-06-0.867E-07 
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RIGHT HAND ROTATION 



LEFT HAND ROTATION 


Source coordinates: in flight system; r^,4>o,x^ in propeller system 
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Figure 5. 



AXIAL DISTANCE FROM PROP-FAN PLANE NORMALIZED BY DIAMETER (X/D) 




Computed Prop-Fan directivity at 0.80 Mach number and 800 fps tip speed. 
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Effect of source placement on computed noise level. Prop-Fan at 0.80 Mach 
number cruise. 
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Figure 7. Unsteady thickness effects for 4-blade propeller at 0.4 Mach number in plane. 
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Figure 8. Contribution of radial loading at BPF for a Prop-Fan at cruise. 
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Figure 10. Effects of unsteady loading terms and non-axial inflow for a 4-bladed propeller 
at 0.20 Mach number. 
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Figure 11. Effect of unsteady loading and angular inflow terms at the peak noise location for 
a 4-blade propeller. 
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PLANE PERPENDICULAR 
TO FLIGHT DIRECTION 



Figure 12. Circumferential directivity of a 0.8 Mach number Prop-Fan in 2 planes passing 
through the hub. 
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270 ° 



Figure 13. Circumferential directivity at BPF predicted with near-field and far-field formulas. 

Prop-Fan at 0.8 Mach number with shaft at 10® angle of attack. Plane through 
hub and perpendicular to the flight path. 


- 86 - 







270 ° 



Figure 14. Circumferential directivity at BPF of a 4-bladed propeller at 0.50 Mach number 
cruise. 5° angle-of-attack, 0.27 diameter tip clearance. Plane through propeller 
hub and perpendicular to flight path. 
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Figure 16. Measured and predicted axial directivity at +7.3“ angle-of-attack. 0.15 Mach 
number, 7=0.717, 183 / 4 = 23 . 7 °). 
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Trends with angle-of-attack for 4 harmonics with microphone at 90°. 0. 15 Mach 
number, 7=0.628, 183/4 = 19 . 9 °. 
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SOUND PRESSURE LEVEL, dB SOUND PRESSURE LEVEL, 



LOAD DISTRIBUTION 




Figure 19. Euler versus UAAP spanwise load distribution for SR7A at 0.20 Mach number, 
7=0.886, Cp=0.87). 
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SOUND PRESSURE LEVEL. dB 



Figure 20. Theoretical and experimental directivities at BPF for SR7A Prop-Fan. 0.20 Mach 
number, 7=0.886, C_=0.87. 
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Figure 21. Predicted and measured trends with inflow angle. BPF for SR7A under flight 
path. 0.20 Mach number, 7=0.886, C^=0.87. 
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Figure 22. Comparison of measured and predicted trends versus inflow angle of BPF for 
SR7A at 0.20 Mach number, 7=0.886, C^=0.87). 
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igure 23. Comparison of the measur^ and predicted BPF directivities at 3 inflow angles 
(PTA boom microphones). 
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Figure 24. Comparison of predicted and measured trends at BPF (PTA boom microphone 
positions). 
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Figure 25. Comparison of measured and predicted trends at 2xBPF (PTA boom microphone 
positions). 
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Figure 26. Comparison of measured and predicted trends at 3xBPF (PTA boom microphone 
positions). 
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Figure 27. Comparison of measured and predicted noise trends with inflow angle at the PTA 
design cruise condition. Unsteady loading based lifting line aerodynamics with 
non-uniform flow-field information provided by Lockheed. 
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ANGLE AT WHICH MIDCHORD LINE INTERSECTS = 
CONSTANT CYLINDER IS GIVEN BY 


tan<^ 
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Figure 32. Sketch for determining angle to mid-chord line at source radius. 
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Figure 33. 
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Sketch for locating source panels for a chordal surface representation. 
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OBSERVER 


Figure 34. Visual and retarded observer coordinates for derivation of relationships between 
systems. 
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